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Von W. Tollmien, z.2. in Farnborough. \ 


Die Untersuchung der Stabilität einer Laminarströmung nach der Methode der kleinen Schwingungen ER, 


% > für große Reynoldssche Zahlen integriert werden soll. In Anlehnung un gebräuchliche Ansätze werden vier, 
0.0... Näherungsintegrale dieser Differentialgleichung definiert und ihre Fehler sowie die, ihrer Ableitungen asymp- 


05 totisch abgeschätzt. Einige Bemerkungen zu dem Randwertproblem und zur Störungsdifferentialgleichung, 


eh, je falls die Grundströmung eine Grenzschichtströmung ist, werden. angeschlossen. 
 rential equation are defined and their errors as well as those of their derivations are asymptotically estimated. 


„Im  coneluding some observations are added concerning the boundary problem and the perturbation 
Bunker differential equation in case that the base current is a boundary layer current. Tr r 
0 Del’examen de stabilite d’un courant lamineux selon la methode des oscillations petites, resulte une equa- 
2. tiondifferentielle,ordinaire et non-adjointe en elle-meme, de la4e ordre, qui est. a integrer assymptotiquement 
Wa quant aux nombres grands de Reynold. Prenant en consideration des disposition d’operation usuelles, 4 inte- _ 
2 grales approximatives de cette &quation sont dejinies et leurs defauts de meme que ceux de leurs derivees sont R 
v övalues assymptotiquement. En terminant quelgues observations sont ajoutees, concernant le probleme de 
la valeur de limite et l’&quation differentielle de perturbation en cas que le courant de base soit un courant 
de couche limite., KEN Me n. 7 
f x i $ Y £ z I } y 
Mccsenogaune YCTOÜYHBOCTH NAMMHAPHOTO TeyeHnnd IIO MeTony Hedolsımmx Kortebaunä 
IPHBOAUT K OÖBIKHOBEHHOMy AUbpepeHIMaBHOMy YPABHEHHIO YETBePTOro paspsına, KOTOpoe 
AON:KHO ÖbIT UHTELPHPOBAHO A3UMITOTHYECKH MIA ÖoANBIIHxX yucest PaliHousna. 
Ha ocHopanmn OÖBIYHOTO MONXONA OMpeneismmorch yerspe IIPHÖHNKEHHLIE HHTErPAJIHU 
9TOrO AuhhepeHnmansHoro YpaBueHun Mu YCTAHABIINBAITCH ABUMIITOTHYECKU UX OIIUOKN, 
3 TaKike ONIHÖRH HX OTBETBJICHHÄ. 
K 3T0My mpucoenuHsIoTcH HeKOTOPbIe 3AMEYAHHA IIO TIOBONY poÖ.temßt IIO ONeHke KpaeB- 
‚(paugeprupoönem) u ambpepenmna,bnoro VPaBHeHNA IIOMeX, B TOM CAYYae, eCJIH OCHOBHOE 
TeyeHHEe ABIAETCH TeYeHHeM TPAHHYYINHX INIACTOB, 


-Die Entstehung der Turbulenz wird auf eine Instabilität der Laminarströmung zurück- 
geführt. Zur theoretischen Untersuchung der Stabilität der Laminarströmung hat man meist 
die Methode der kleinen Schwingungen benutzt. Reynoldssche Zahlen für die Stabilitäts- 
- grenze wurden nach diesem Ansatz 1929 von mir. für die Strömung längs einer tangentiell ange- 
4 _ -strömten Platte 1) berechnet. An diese Arbeit schlossen sich zahlreiche und ausgedehnte Unter- 

suchungen von H: Schlichting und J. Pretsch.an. Es sei etwa erinnert an die Unter- 

% suchung von H. Schlichting über das Verhalten von Strömungen mit Wärmeschichtung ?); 
= die Ergebnisse dieser Rechnungen stehen in ausgezeichneter Übereinstimmung mit genauen Ex- 
'perimenten von H. Reichardt. Von den Untersuchungen von Pretsch sei beispielsweise 
"angeführt seine Übertragung der Stabilitätsdiskussion von ebenen Strömungen auf rotationssym- 
= metrische 3). Dadurch wurde ein wertvoller theoretischer Beitrag zum klassischen Stabilitäts- 
; problem der Hagen-Poiseuille- Strömung im Kreisrohr geliefert. Diese Erwähnungen 
greifen nicht mehr als Stichproben aus der Fülle des neuartigen Stoffes heraus, der durch Schlich- 
ting und Pretsch erschlossen ist. . | 
‚Alle diese Stabilitätsrechnungen sollen asymptotisch für große Reynoldssche Zahlen 
gelten. Einige andeutende Bemerkungen darüber, in welchem Sinne die Störungsdifferentiäl- 


We Woltmien; Über die Entstehung der Turbulenz. 1. Mitteilung, Nachr. d. Ges. d. Wiss. zu Göttingen, 
Math.-phys. Klasse 1929, S. 21 bis 44. Die jetzt vorgelegte Arbeit kann man, wenn man will, als verspätete 2. Mitteilung 
"zu dieser Untersuchung auffassen. 

3) H.Schlichting: Turbulenz bei Wärmeschichtung. Z. angew. Math. Mech. Bd. 15 (1935), S. 313 bis 338. 
- 3). J..Pretsch: Über die Stabilität einer Laminarströmung in einem geraden Rohr mit kreisförmigem Querschnitt. 
Z. angew. Math. Mech. Bd. 21 (1941), S..201 bis 217. 
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ntialgleichung ebener laminarer Strömungen bei hohen 


. führt auf eine gewöhnliche, wicht selbstadjungierte Differentialgleichung vierter Ordnung, die asymptotisch _ : 


\ x - r RER a ee Py- p D EN h RAN RE 2.) 5 h \ e® 
The examination of stabiliy of a laminar current \referring to ‚the method, of small vibrations is re- 
-  sulting in an ordinary not self-adjungated differential equation of ihe 4" order, which is to be integrated. for ı 
00. high Reynold’s numbers. In connexion; with current arrangements four approximate integrals of this diff- 
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N ‚asymptotischen Integrale. ae Bu 


1, f Die , Störungsdifferentialgleichung. "Die ebene‘ laminare Strömung Bbe nur eine ei - a 


es schwindigkeitskomponente U in der x-Richtung, die nur von der Koordinate y senkrecht zu = 
abhängen möge. Die Stromfunktion einer störenden Partialschwingung sei mit p(y)eirt@—en 


5 . angesetzt. a ist bis auf den Faktor 2x die reziproke Wellenlänge, c ist die Präsens 2 
der Störung, wenn wir. crein reell annehmen. Dies trifft zu für die neutralen. Schwingüngen ander 
‚ Stabilitälsgrenze, für die wir uns interessieren wollen, Die Differentialgleichung für die komplex-. R 


wertige Amplitude p wird aus den Navier-Stokesschen Differentialgleichungen unter. Be- 


SO RABN. gu lineare. Glieder en 2); und lautet in dimensionsloser. ei ns BEER N 3 
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Die linke Seite dieser Differentialeläan entspringt aus den Trägheits- und Druckatiedene der ur 


Bewegungsgleichung, die rechte Seite aus den Reibungsgliedern. Striche bedeuten stets Ablei- 
‘ tungen nach. Y.. RisteineReynoldssche Zahl, die mittels einer 'charakteristischen Geschwin- 

“ digkeit und einer charakteristischen. Breitenerstreckung des Geschwindigkeitsprofils gebildet ist. 
‚Diese Differentialgleiehung soll sogleich etwas vereinfacht werden, indem von den Höäbumgsficdernd 


. auf der rechten Seite nur das von der höchsten Differentiationsordnung mitgenommen wird. te 


kann nämlich’statt Gl. (1) nach leichter Umtortns ran 


BR a. ’r 2 (v" 22 =) a RAcH ehe - | 2 
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Nun soll die Untersuchung dieser Differentialgleichung asymptotisch für große Reynoldssche 
Zahlen, oder, worauf es hier ankommt, für große Werte des zusammengesetzten Parameters « R 
durchgeführt werden. Da bedeutet es nur eine winzige Abänderung der Phasengeschwindigkeit ec 
bzw. der Funktion U”, wenn wir statt A a) die schon VEREDIGIE vereinfachte ar 
gleichung 


(U —)(d’ —a29) -Ug = —— ME 


nehmen. Übrigens Äh ‘de die asymtotische Diskussion der vollständigen Differentialgleichung (1) 
keinerlei zusätzliche Schwierigkeiten, sondern nur etwas mehr Schreibarbeit bedingen. Es würden 


eben Glieder mitgeführt, die bei dem schließlich gewählten Annäherungsgrad. nieht ins Gewicht 


fallen. 

Für das Arbeiten a, der Differentialgleichung (2) spricht noch ein anderer Umstand. Die 
Annahme, daß U nur von y abhängt, trifft nämlich fast nie genau zu. U hängt meist auch von x: 
ab, und die laminare Grundströmung hat demnach noch eine Komponente V in der y-Richtung. 
Schon 1929 (a. a 0. 5.40) habe ich darauf hingewiesen, daß bei einer sehr wichtigen Klasse von 
Geschwindigkeitsprofilen, nämlich bei den Grenzschichtprofilen, der Einfluß von & auf U nur. 


scheinbare Reibungsglieder niedrigerer als 4. Differentiationsordnung in die Störungsdifferential- 


gleichung einführt. Nähere Ausführungen über diesen Punkt wird $7 bringen. Da diese schein- 
baren Reibungsglieder ohnehin gestrichen werden, ist es zweckmäßig, von vornherein nur das 
höchste Reibungsglied zu berücksichtigen, wie es in Gl. (2) geschehen ist. 


1,2 Einführung der asymptotischen Integrale. Das mathematische Problem, das wir uns 
stellen, ist also die asymptotische Integration der Differentialgleichung ( (2) bei großem « R, Eine 
wesentliche. Verwicklung kommt dadurch, hinein, däß man eine Nullstelle von U —c im Innern 


*), Für den Leser, dem diese Stör BER RER nicht geläufig ist, sei auf $ 7 verwiesen, wo eine Herleitung 
unter sehr allgemeinen Voraussetzungen gegeben wird. 


= SUR 1. Die Störungsdiferentiaigleichung Einführung ihrer “ Y Br 


_ gen ganz auf derartige Grenzschicht- 


- örterungen für den einfacheren Fall 


des I auf die Randbedingungen (vgl. $6) 


d etes mit Rücksicht 
genüg 


y= © beibehalten werden soll. Eine genauere Berücksichtigung des asymptotischen Überganges 
von Uin U„ dürfte man nicht als lohnend ansehen, wenn man bedenkt, welche sonstigen Ver- 


‚ nachlässigungen bei der AufstellunG 


der Störungsdifferentialgleichung 'ge-. 
macht wurden. Ber en 
Wir werden unsere Betrachtun- 


profile einstellen. Besondere Er- 


symmetrischer Kanalprotile sind dann 


_ kaum noch notwendig. In Bildlist #7... .0 a2 13 0 23 0 


R] 


das Grenzschichtprofil aufgetragen, Bild 1. Geschwindigkeitsprofil an der enter angeströmten 


. das sich bei der tangentiellen Anströ- Platte. c = 0,1328 bzw. 0,3938. 


mung einer ebenen Platte ergibt.  _ 

Das erste Paar asymptotischer Lösungen der Differentialgleichung (2) könnte man so auf- 
zustellen versuchen, daß man die mit dem sehrkleinen Faktor, multiplizierten Glieder von Gl. (2) 
überhaupt wegläßt, also mit der reibungslosen Störungsdifferentialgleichung 


EN. BR a Er en re 
- rechnet. Da 
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g® der kritischen Stelle U = c einen Pol erster Ordnung hat, so liegt dort 


eine singuläre Stelle der Bestimmtheit vor. Ein Lösungspaar von Gl. (3) drückt sich mit %, und 


SP, als zwei regulären Funktionen so aus: 


ra | anzunehmen hat. Es 
2 gt aber, nur eine einzige solche Nullstelle in dem Gebiet vorauszusetzen. Von ihr aus wollen 
wir die Querkoordinate y messen. Größen an dieser kritischen Stelle sollen mit dem 
© Index 0 gekennzeichnet werden. Die positive y-Richtung soll von der Wand fort, nach außen 
hin weisen, und es soll U, >0 sein. Ferner soll U‘ nicht verschwinden; U, = 0 wäre ein Aus- 
-  nahmefall, dessen Erörterung wesentlich einfacher wäre. Das Integrationsgebiet beginnt an der 
Wand, wo y gleich der negativen Größe y,„istund U verschwindet. Beisymmetrischen Strömungen _ 
durch Kanäle erstreckt sich das Integrationsgebiet bis zu einem positiven Ya, wo U sein Maxi- 
mum U„ annimmt, also U’= Oist. Bei den besonders wichtigen Grenzschichtprofilen reicht das 
Integrationsgebiet an sich bis y— ©. Die Grenzschichtprofile seien nun insofern etwas idealisiert, 
_ alsder Maximalwert Um von U nicht erst asymptotisch bei y— oo, sondern schon bei einem end-. 
‚lichen Wandabstand, der einen endlichen positiven ya entspricht, erreicht und dann konstant bis‘ 
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Dabei beginnen die Potenzreihen für P, und ®, so: 
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Für positive yistin Gl. (5) log yrein reellzunehmen, für negative y dagegen gleich log |y| — iz, 
wie von mir 1929 (a.a. O. S. 29) nachgewiesen wurde und später nochmals erhärtet werden wird. 

Während o, nach Gl. (4) auch im kritischen Punkt regulär bleibt, wird o,, also die «-Kompo- 
nente der Störung, dort im Widerspruch zur Methode der kleinen Schwingungen logarithmisch 
unendlich, erst recht wachsen die höheren Ableitungen von 9, dort über alle Grenzen. In einem 
schmalen Gebiet um U=c, d.h. um y=0 muß daher jedenfalls die Reibung berücksichtigt 
werden. Setzt man in diesemschmalen Gebiet o,annähernd gleich 1 [vgl. Gl. (5) und (7) bei y=Q], 
U—c gleich U, y, U’ gleich U,, so erhält man für die an Stelle von 9, infolge der Reibungs- 
korrektur auftretende Lösung g(!) näherungsweise die Differentialgleichung: 


‚ [22 7 % 111 : 
U, yo» er ELLE RN MORE .(8). 
Es sind g(V solche Bedingungen aufzuerlegen, daß es außerhalb eines schmalen Gebietes um y = 0 
in die reibungslose Lösung 9, übergeht, was in der Tat möglich ist. Weiter wollen wir 9, und go 


zu einer einzigen Funktion p, zusammenfassen, die für das ganze Integrationsgebiet als Näherungs- 


ar, 


rang 7 = a nd BEZ are BE a ad sin 
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lösung gelten soll, also für hr eine |y| mit 90 zusammenfällt, während sie für größere Is sich 
9, nähert. Wie in’ 4,1 gezeigt ala wird, gelingt dies durch den Ansatzı 


BB) 2.2... BEN; 


Außer den bislang aufgeführten a Integralen g, und 9} benötigen wirnatürlich 
noch zwei weitere. Betrachten wir als einfaches Modell von Gl. (2) die Differentialgleichung 


v0 


Yo" —ioaR@"—p)=0, 


sosieht man, daß nntohisch auch zwei Integrale auftreten, die bei zunehmendem a R unbegrenzt 
wachsende Differentialquotienten aufweisen. Daher ist für dieses Paar von BSNTEBA TE LIE Inte- 


gralen die Differentialgleichung & 


ga" —ioRg’ = 0 
maßgebend. Entsprechend setzen wir für das zweite Paar rn 0 Lösungen 3 und OR 


von Gl. (2) die Differentialgleichung 
ON NER NR 0 ee are (10) 


an, indem wir nur die höchsten Differentialquotienten aus Gl. (2) berücksichtigen. Freilich hat 


Gl. (10) auch noch eine beliebige lineare Funktion als triviale Lösung, die zu unserer Leit- - 


anschauung nicht passen würde und daher durch geeignete Anfangsbedingungen auszuschließenist. 

In Anlehnung an frühere Vorschläge wäre man versucht, das erste Lösungspaar ,, p% als 
„langsames‘‘, das zweite Lösungspaar o,, p, als „Tasches“ zu bezeichnen. Eine solche Benennung 
wäre indessen wohl doch nicht sehr glücklich, da } sich in der Umgebung von y = 0 auch rasch 
verändert. 

Unsere Aufgabe können wir jetzt genauer so fassen, daß die Abweichungen der asympto- 
tischen Näherungen 9,, 9), 9 und 9, von den zugeordneten exakten Integralen 91, 11, Yım und Qıv 
abzuschätzen sind. Die Asymmetrie in der Bezeichnung der asymptotischen Integrale können 
wir getrost in Kauf nehmen, da go} in der Tat eine Sonderstellung inne hat. Allgemein sollen 
römische Ziffern als Indizes der exakten Lösungen dienen, während arabische Ziffern zu den 


Näherungen gesetzt werden. Ein zugesetzter Index m soll das Maximum u betreffenden Größe 


bedeuten. 

1,3 Exakte Lösungen der Störungsdifferentialgleichung für y > y,. Zunächst seien einige 
exakte Lösungen angegeben, die in Spezialfällen sofort hingeschrieben werden können und uns 
bei der Fehlerabschätzung nützen werden. 

Bei Grenzschichtprofilen vereinfacht sich die Störungsdifferentialgleichung im Gebiet 
ana. W0.0U— D.ist, zu 

en pP" —iaR(Um—) (RO) —=0. 
Mit dem Ansatz 9 = eV bekommt man für die möglichen vier Konstanten % die Beziehung 


lutra 44.02 
= 1 1+ — ——— |... .2 202. : 
fat ge VERER (11) 
Daraus bekommt man a für das erste Konstantenpaar 
Ve ERROR 
Mara, a SE OR (12). 
Man darf also zwei Integrale 
De 
Der, 92 = — Eimay Re em Abe EEE RER) 
9) 


setzen, solange yS* x Ist, wohlgemerkt nicht bis zu beliebig großen y. Für das 


zweite Konstantenpaar ee man angenähert 


bus F ViaRUm—g| ee 
a F io R(Um—c) an. n RA ERTL (14), 
und man kann entsprechend zwei Integrale i 
1+i 
P; 3 eXp et Rn) SE es RER che (15), 
l+i. 
Ps exp | 2 Y& R(Um—c) (vH) Bu 20 Re Rn 0 A (16) 
setzen, solange y< VaR(Um—o) ist. 
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‚eine neue abhängige Variable F eingeführt: 
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Dr: Asymptotische Diskussion zweier Hilfsfunktionen fı und fır. 


2,1 Einführung neuer Veränderlicher in die Differentialgleiehung der Hilfsfunktionen. Wir 
wollen ‘zunächst Hilfsfunktionen f bei großen & R studieren, die der Differentialgleichung 


Te BR EN ER Be NEE ER de et) 


genügen. Diese entsteht aus der Differentialgleichung 1 G]. (10)?), wenn man ” proportional zu f 
setzt. Der Vorzeichenwechsel von U —.c erlaubt esnicht, fertige asymptotische Lösungsdar- 
stellungen Zu übernehmen. Wir werden die asymptotische Diskussion von Gl. (1) auf dieLösung 
einer Differentialgleichung von. der leicht i 
zugänglichen Form 


hinausspielen. Zu dem Zweck wird statt y 
eine neue unabhängige Variable :Y, statt f 


3 2/3 Den 
FI) m WILD eaNe, 2 
() {UV f 
rs RR IR AR 8), 


Diese Transformation ist von dem Para- 
meter & R unabhängig. Ferner ist 


dY U—c 
For N) 


0 


dy 
SER, 3 Bild 2. Y und —— für das Geschwindigkeitsprofil an der 
stets positiv und in der Umgebung von dy 


y— O ist tangential angeströmten Platte bei ce = 0,1328." 

Er 2 2 & 

2 ZUR) 
173028 

Y ist danach eine monotone Funktion von y, die in der Nähe von y= 0 in yübergeht. Das Inter- 

vall'von Y erstreckt sich von Y,, bis oo entsprechend dem Intervall von y, das von %, bis oo reicht. 


10% 
ya BR ENG 2 


da), 
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Bild3. Y und E für das Geschwindigkeitsprofil an der 


Y 
tangential angeströmten Platte bei ce = 0,3938. 


2 


RS d} 
Die Bilder 2 und 3 zeigen den Verlauf von Y (y) und dr für das Geschwindigkeitsprofil an der 


tangentiell angeströmten ebenen Platte (s. Bild 1). In dem Bild 2 ist c = 0,1328, in dem Bild 3 
ist c= 0,3938, wobei U„ = 1 gesetzt ist. Für y > y, ist 
5) Auf Gleichungen aus demselben Paragraphen wird durch einfache Angabe der Gleichungsnummer verwiesen. 


Bei-Bezugnahme auf Gleichungen aus früheren Paragraphen wird die Paragraphennummer noch vor dieGleichungssnummer 
gesetzt. 


en 
bleibt stets beschränkt, Her y- ya 0 wird. | 
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wie man am Rn. un GI. (5) erkennt. nr N geht K 2. y*' oder mit Ynach 0. B 
Danach liegt es nahe, Re ag IategsnlE = und a von g: (7) varuetır aus re 
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ar —ia «RU, Yr,, BZ ni 3 He 2 ah .®. 
Be 3 zu berechnen. Mit! 0% ee... 
es Ne Re RE NE = (aRU,)j8 ee RR ee 
ae nen 

erhält man aus Gl. (9) NEE | FE ESTER 


aan 20 N RN ee Er 


H erstreckt sich von H,, bis oo. Die H an k e REN erster und zweiter Art von deu ER 3”. R 
Ordnung 4 [3 gestatten, die ea yon Gl. (12) aufzustellen: N | I 


F,= Hu2 u) I ap]. ee un K a 
r,— ie em] BER Par . am) 


Um bequeme Formeln für dieHan n kelschen Funktionen bei großen |H| zu bekommen, re GL 
F, und F, in der unteren H-Halbebene betrachtet. Für positiv reelle H ist (#. )#'2 also gleich Fee 
eöltin 912, ie negativ reelle H gleich a |2 mit positiv reellem 5 H 32, Für kleine | A| wird- 20 0 


m re R 
’ sin 1/3 et 3 3 

) isin | 5 (5) sin I. 

j in } "Min 
313. 2% , En 
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IT bezeichnet die Fakultät. Nach einem elementaren Satz aus der Theorie der Differentialglei- 
chungen ist die Wronskische Determinante des Fundamentalsystems F, F, von Gl. (9) 
konstant, und zwar ist nach Gl. (14) 
» dr, dF, 6 | Be insp 
j m ehe bi RG EN Se dealer ee BEERle LE ENE e 
AH aH! FR | z RR (18). 
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Nach den Formeln für dieHankelschen Funktionen bei großen IH (vgl: etwa G.N. Watson: “ 
A treatise on the theory of Bessel Funktions 7. 2 p..196 sq., Cambridge 1922) wird: 


> Fım K e = meutexp| 1, arme] - ; R ' 

2 | N RR RT 16a), re 
| 5 385 ( 

Ä RT EA N al ee) (H-212 ‚ 

i | Erehe ı)H 1608 °# rO(H } hr 
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Oistdas Bachmann-Landausche Symbol für asymptotische Größenordnung. F, klingt 

danach für wachsende positiv reelle H exponentiell ab, für negativ wachsende reelle 4 nimmt 

es exponentiell zu; denn für negativ reelle H ist 432 gleich ö| H]??. F, verhält sich gerade um- 1 
gekehrt. | 


2,3 Abschätzung der Fehler von F', und F',. Es gilt nun festzustellen, inwieweit F, und F, 
als Lösung von Gl. (9) gewisse zugeordnete Integrale F, und F}ı der exakten Differentialgleichung (7) 
asymptotisch darstellen. Für den Fehler #,—F, bekommt man nach ‚Gl: (7), (8) und (9) die 

. Differentialgleichung m 


ee 
ms meu MA r)= En 

oder | 
drlR a, “ | 
a) Be EN 1 1 RE (17). 


Mit einem auf J. Liouville zurückgehenden, öfter angewandten Kunstgriff fassen wir Gl. (17) 
formal als inhomogene Differentialgleichung für #; —Ff', auf, deren homogener Teil durch dasselbe 
Fundamentalsystem wie.Gl. (12) befriedigt wird. Mit Rücksicht auf Gl. (15). erhält man: 


Mer MEIRIMAH NT RE NAHEN Fer (18). 


Die Integrale sind längs der reellen Achse zu erstrecken; die Wahl ihrer unteren Grenzen erfolgt 

später, sieist gleichbedeutend mit der Festlegung der exakten Lösung Fr, mit.der die Näherung F, 

verglichen werden soll. Würde man beispielsweise,beide Integralgrenzen nach demselben Punkt 4, 

- legen, so hätte man über die Anfangsbedingungen für Fin dem Sinne verfügt, daß Fj bei H, 
mit F, übereinstimmt, desgleichen F, mit F\. 

Für F, und F, führen wir eine andere Schreibweise ein, die an deren Form Gl. (16) für 


große |Z| anknüpft. Das Symbol o, (H) möge bedeuten: 
u au IS N (19). 
o,(H)=|Aj" für |24>1 J 
Dann erhält man die Darstellungen 
Po.) exp| + us DE ern Na et (20a), 
2 | 
M=0_14(H) exp | un (1-+:) uw RE RE RL: (20b), 


‚wo 9, und p, bei beliebigem Wachstum von H beschränkte Funktionen sind. Entsprechend setzen 
wir die exakten Integrale an: 


"3 
FL =0_1,(H) exp Ei „(1+:) Han] DELETE ee Pe REES GR 


ee 
2 ; 
Par =6: 1/4 (H) exp | ; r (1 + i) a yanı (H) EEK TR (21 b). 


BEN zunächst N mn a wir aber als untere Grenze des zweiten Inte x 


Für 7 die ee der Beträge a Bee eis von ‚Cl. Le EN die 


‘das Maximum von H, nämlic 0, SO} 
trag der Exponentialfunktion des Integra 
oberen Integralgrenze H seinen größten Wert 

‚ denn dieser Betrag ist für positiv reelle wie an 


exp (—2 „12 se) und für negativ reelle Kl gleich ,_ 


zweiten Integrals in Gl. (22) Fa man daher 


vor .das Integralzeichen 


\oxp |—2 „12 a; ru ‘) ma 


Ve 


Grenze des H-Intervalls, zusammenfallen. Nach ‚dieser Festlegung, aus Gl. Er und vn 


; FE 


N en Hufen, ann] Bumam >) 


Ar, H 
\ { FEB 
Pr 591 = in Ironman : 
v . (24). 
— Pzexp Ba (1 + ya | Ro-men Be 21. (1 +) ze |pı nd) 


\ Bessiehnst man mit dem Index m die in der betreifenden Beträge, : so liefert eine rohe Ab- 


schätzung von Gl. (24): 


} |?rlm &IPpılm + edle { Tonikjan + op Klam). 
6 ” Ha 


Y 
Indem ja statt H eingeführt wird, bekommt man 


’ S oo 
| ENT 1 
N lvrlm<|Ppılm or 5 ?ılm|2ılmIPelm) 0-12] K]d y 
also | “ 
» f 
In ee Alm Ne ee 5). 
1 animal] BERBIEIG E | 


Das Integral Auf der rechten Seite dieser Ehaleichnut konvergiert zufolge des ee Ver- \ 


haltens von K bei großen Y. & muß natürlich so klein sein, daß der Nenner auf der rechten Seite 
von Gl]. (25) positiv ist. Benutzt man diese rohe Eingrenzung von ?1, um die rechte Seite von 
Gl. (24) nochmals abzuschätzen, so bekommt man als Serfeinerie Fehlerschranke: 
Ar Inka) o-unlKld r 
Meise un DE ee Te (26). 
ganzen ma, d Y 


als absolute‘ Schranke der sag a a FR 


3 en Er) || ziehen und iR dem Betrage der schon ohnehin vor 
dem Integralzeichen stehenden Exponentialfunktion 

Br RE zu 1 zusammensetzen. Die noch offene untere Grenze 

- des By Integrals kann’ man Heinen, beliebig wählen. Wir,lassen sie mit A,, der unteren 


ep[? Ei (S. Bild 4). Bei der Abschätzung RE = Br: 


er 


1. = 
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Es genügt für unsere Zwecke, o_ 1a in dieser Formel durch sein Maximum 1 zu ersetzen; wenn 
man auch durch genauere Berücksichtigung des Verlaufes von o_1)g den Fehler von p, von der 
Größenordnung, e auf die Größenordnung 3/2 hinunterdrücken könnte. 


Als Bestimmungsgleichung für den Fehler von F, legen wir nach entsprechenden Über- i 


legungen 
/ Aer H H 
Fun, {Ps SER MdH-mMlKr nom. 227) 
Bw Se RR | 
zugrunde. Man erhält daraus als rohe Schranke: | 
rim [Palm En Be Ente (28) 


Een 
1 |Pılm|Pe m Je-wlRlar | 
‘und als verfeinerte Fehlerabschätzung: | Br 
S oo 
2 |»1lm|P2] m/ o-1e| Klar 
Im— |< BEN PEN re (29). 
\ EN 
I IPılm| Palm [ 0.12) K|AY 


2%,4 Wronskische Determinante ‘von fr und fir. Für spätere Verwendung wird die 
Wronskische Determinante ; 


a. A 
dy’  dy 
benötigt. Nach Gl. (3) ist ’ 
dfı EN RR dFır 
dy Tag A=-7Gy II ray F} Hi BT N A Le Une: (30). 


Wegen der Gestalt der Differentialgleichung (7) ist diese Wronskische Determinante für alle Y 
konstant, so daß es genügt, sie für ein Y, etwa Y,,, zu berechnen. Nach Gl. (23) und (27) bekommt 
man: 


d# ET 
ay #: en )lit | er, mar). 
Yy 


dF; dFı Sn 
Bee are 


Nach GI. (15), (20) und (21) wird 


dF dr 6% ef R 
tn m le 0-12 KpıpudY) RENNER ER ARE (31). 
Y 


ET 


Mittels Gl. (28) wird das zweite Glied in der Klammer abgeschätzt. Danach gilt die Näherung 
dfFı dFır 65 
dr Fr — dr Fj Lem En Ne re ee ee Alain nei idee 


bis auf einen Fehler, der höchstens vom Betrage 


[0,0] 
I2lmlvalm/ 0-12 Kar | 
Yıv Re AR Kr a ie Pauls Mist hrs (33) 


EI 3 
1 Be l2ılml22lmf 0—1/2 |KjaY 
! Yu 


ist. Wenn man hierin o_17» nicht durch die obere Schranke 1 ersetzt, sondern den. definitions- 
gemäßen Verlauf berücksichtigt, so kann man einsehen, daß dieser Fehler sogar mitye klein wird. 


- dF dF 
2,5 Fehler der ersten Abteilungen von F, und F,. Sodann soll der Fehler von AH und TH 
abgeschätzt werden. Nach Gl. (23) ist 
dM—F). 2rjdef dr, t \ 
Bau lm man Ss elxmrdaH.. %... 34). 
er en ee ze En 
oo 


tolgt.. Habe ist die N so u daß cine ochmalige 
‚a. Ne) liefert. ie ae rt: Br $) ER ar £ 


en: | . a "a. Eu a4) ver +0 u: “ a Ki 
Ei: Also man En ER: Er Sr an # EB EN = ; 
ER: 2 au.) exp Burn] Pa. 

Bi BD & dR, 


PR Re N “ ER Er a an (A) exp Ba om] Pac) SEE ; 


setzen, wo Bun und‘ Pa ‚ beschränkte Funktionen sind. Entsprechend wird 


NT: dr; — ua ( (2 exp ns 143) mil et Sa 


dB. 
am, 3 
en lee une er Ba pet RN 
angesct Damit wird aus GI: (34) 
| 2 ve 52 
ir ArD— 1 lau Be 3 
h i w.' % ; $ } ‘ 2 
H z 2 Sr 
\ Bi 2y2 SEHE RN 
| a4] [urn a+ 9 mel nn.az 
Also wird EEE | READER, Sa 


ö Yr - \ \ © b ; 7 , 
Ir rn |< tlPumloılm| Pelm) 022 |K|AY+ |p.ı Im Pl [ulm | al yH Eee 
Darin ist noch die obere So von |Pı "a nach Gl. (25) EUEULaE SI Entsprechend wird 


|Pırı — Pa |< T ou lm! Dil Pal oe XIAY Hl Bull o-ı2|Kldr} (40), £ 


4 \ 


wobei |pır|m nach Gl. (28) in die rechte Seite einzusetzen ist. 


2,6 Spezielle Näherungsausdrücke für fi und fir. Es mag angenehm sein, gewisse Näherungen. A 
. für die Lösungen f(y) der ursprünglichen Ditferentialgleichung (7), auch einmal unmittelbar vor 
Augen zu haben, ohne den Umweg über F(y) gehen zu müssen. Aus Gl. (14) folgen bei y=0 


die Näherungen: . 
m \ im 
ae use‘ 
fe ee er 
11(—,) sin x (4 )isi 3%: 
| 3 EN BEE 


| u i 
Die Umrechnung von Gl. (16) liefert für solche y, bei denen YaR j [ Y U—cd y| groß ist: 
0 


.der Zweig e a7 (e— U)-Vt zu nehmen, wie eine nähere te Erörterung ( der Herkunft dieser Wurzeln. 


aus den Formeln für F(Y) AerEh \ 


3. Veh ehschschitsung von Y5 und Du 


x 3,1 Detikition von 9 und 9;. Verhalten von 9, bei y>Yy,. Gemäß 1 Gl. (10) ind 2 Gl. “ 
Ben wir o, ra zu Ir p, proportional zu fi an, und zwar wollen wir 


UNI NE LET ENTER NER 
E n P; Fiah: N el: an ene h I  RRRERNER $ wie Ma (1) 
' nehmen, so daß mit BLEI \ E 
| F} \ IR 
4 N = PR EEE Er (2) 
/ d? -1,.d? 
in = /fı 2 gilt 


wird. Über den Fehler, den wir bei Ersatz von fi durch f, und von fır durch fs begehen, wissen wir 
„aus dem es Paragraphen 2 Bescheid. Vorbehaltlich 'späterer BegrimlLung setzen wir weiter 


n 
1 1 
P, = ZZ = [kan u Art FE Mei dt cal any Bude. Herz alıkhe re Zee (3a), 
H 2 nn 
®3 a) Snavas= | Snanan a re ic er (3b), 
[m N NEN SEE er f . (da), 
Yw My 
Yu a Tank RD EN (#b). 
‘ Yw Yw um) N 


Dabei haben wir uns auf Grenzschichtprofile beschränkt, bei anderen Profilen könnte 'man die 
Wahl der unteren Integralgrenzen in Gl.(3a) und (3b) anders vornehmen. 

Die Integralgrenzen i in Gl. (3) sind so eingerichtet, daß p,im Gebiete des konstanten U= U„ 
mit guter Annäherung i in die dort definierte en Lösung 9,; [vgl. 1, Gl. (15)] übergeht, 


. wie wir jetzt zeigen werden. 
Nach Gl. (1), 2 (3),2 (4) und 2 (16a) u wenn man fr durch f, mit einem Fehler der Größen- 
ordnung ersetzt, für ah H 


19. 
El N u Ar % 2 3/2 | — 3/2 5 
ara, Puls) er 12 ar) Helft 708-3)... 0). 


- Das Streichen des Restgliedes in diesem asymptotischen Ausdruck bringt ersichtlich einen weiteren 


Fehler, der klem von der Größenordnung e ist, wenn 
3 NE ER PETE ER RECENT (6) 


ist. Diese Voraussetzung kann man für y>y, in den praktisch vorkommenden Fällen ohne 


Einschränkung hinnehmen. -Man vergleicht nun diese Lösung mit der Speziallösung @,,, zu der 


eine freie Konstante d, hinzugefügt sei: 


_ Während bei. Bl ln U—c he Wurzeln i in ee Ausdrücken Re zu he sind, ist ” x © 
bei pegaliv reellem U —e für YU=ec U—-c.das negativ imaginäre =. und ‚für i 0: re | 


2, Um Übereinstimmung von Gl. (7c) mit Gl. (5) zu erhalten, muß man 


IR wählen, wenn man einen rn Fehler der RE E Trial) Setzt man A wie wir { 


‚anzubringen ist, die durch genauere Bestimmung von k, nach 1 Gl. (11) zu ermitteln ist. 


wofür nach 2. Gl. (4) 


er u er gi Be u) wi 
Fu. Roi ul i 


„Iti 
IE 


dp = &% ior (Un) exp 2 


Sn . Ann 


er) |. 


rt 
\ ; 


Er E 


ji RR ze 08 (Um —eyBl4 |, fr). ..® (8) B 


es in Gl. Er eigen 


= man. 28 he 


so sieht man aus dem für YZYy gültigen Ausdruck Gl. (5), daß man bis auf Fehler der Größen- 
ordnung e auch Übereinstimmung von 9, mit d, 9, und von 9, mit dy ss ‚bekommt: 


\ 


= dpsll+0(e] für Y= Ya > a EEE SE EIN 


Natürlich verliert diese Aussage diesen Sinn für ganz große %, wo an ,, selbst eine Korrektur 


- Mit dem Vorstehenden ist selbstverständlich noch nicht gezeigt, daß die 9, zugeordnete BR I 
Lösung gr für y> Y, in eine abklingende Lösung übergeht. In 86 nn eh besondere 
Erwägungen. 


3,2 Darstellung von P; und 9. sowie Ihe ersten drei N Zur Vorbereitung der 
Peerabschätzung, leiten wir Darstellungen von 9, und 9, sowie ihren ersten drei REN 
ab. Nach Gl. (1), 2 (3) und 2 (21a) ist 


me are) min. 


| v,y\u Er 
en u) exp Farm] Psa(H) » . 5 (10) 


gesehrieben RER kann. 93, ist beschränkt. Aus Gl. (1), 2 (3), 2 (4), 2 (16a) und 2 (35a) folgt 
für große |H|: - 


; 1/4 
Er site u a 2 En 12a 14 a] OH}. . 1). 
Daraus gewinnen wir die Darstellung 
| et U,y\-ıa SE To | 
u near. 13) 
mit beschränktem 933. 


Um ähnliche Darstellungen für , und , zu gewinnen, benötigen wir eine Integralunglei- 
chung. Es sei w eine Funktion von Y, die stets größer als O ist. Außerdem bleibe auch bei be- 


liebig wachsenden Y, wie sie bei Grenzschichtprofilen auftreten, "die logarithmische Ableitung 
von w beschränkt: 


d log w 
| dY su u Be LE (13). 


DRM 
Für w werden später Potenzen von mr gesetzt. Ferner sei 


ER 


3 | pam R 
; ; far das, age der Iran der beiden in Gl. os ergibt sich: 
Een u = d m U2, 


2 ven man GI. u) an) und as) Br Damit Holst aus En as: 


e- ; „a Rn A 2m) f Hrw exp ( u) dH s namens il am), | = =. 


oda die: Be: Ungleichung für 221 Bestätigt ist. 
| Beim zweiten Sonderfall sei. 


Bi. Dr | . ER | 38 RE DREIER 
E Au Ki J ER H< An > Hn: EL Min =: ve een 5 il a % an 49). 
BE 5 urn, Zeitegung des re erhalten wir SE j TER \ | 

m 2 SE Eu re N 

un ; | Be 

>: Für ie Abschätzung des ersten eh, a der rechten Seite von GI. (20) schen v wirvndr En. 
0 folgenden apa aus: | b 
Be —_|gr-12; | 10 
2 :B, ereor exp > un) d| A| — eh Ka Am AP Aw exp und 
bi Be 3 
F RE NE 
Be | er 


Das a 5 ee des zweiten IS ass zu dem des Stan tal ist 


anf Ze ne it ae 


_. 


. (22), 


wenn man Gl. (13) und a berücksichtigt. Damit erhält man aus Gl. (21): 


L —l2e ") [ Vase (21a) apaıe yauınn-"ewexo (2 pe) 3 3 | 


[#n| 


% Di. Bas Ah Er der. EN Seite von Gl. A 
ersten Sonderfalles als beschränkt. en ‚kann daher mi a | 


n Zieht n man Gl. .(23) and (24) zusammen, s0 hat ER ss a 


* die untere Frtegralkreuze zwischen H, und 1 BE keinerlei Schwierigkeiten, da sowohl das ne er 


a in Gl. (15) als auch nu 


\ | few on Be u+9 an 2. 


a: 


N 


I ar hs 


en a: ie s 2% a 


® / ne 


m el a) 1 


\ 


Damit ist die Ungleichung Ach im uveien Sonderfall bestätigt, 
Schließlich: bietet die ‚Bestätigung in dem Gebiet zwischen den Helen "Sonderfällen 


exp [5 BR de D sm dort endliche, von 0 Vonkchrenehe oh 


- Werte annehmen. . ae | RER 2 f ra 
‚ Für s erhalten wir nach Gl. ) und ( (a): \ % 
j E H j . - x + Y 2 
a a ee 
LER Jay Ir E Jay r > 
. Nach GI. (10) und der eo (15) ne ‚wir also die Darstellung: ER c E “ 
| a1 DPA we 
Ps a 0 - 3/4 e we a NR | Pa(H) . 2.0... (26), 
wobei P3;, eine beschränkte Funktion ist. Entsprechend erhalten wir: e s 
U,y\5% 2 
nee area ....2... en 


mit _beschränktem 73. 


Für p, und seine drei ersten Ableitungen bekommt man mit Hilfe der zu Gl. (15) analogen 
Integralungleichung 


H 2: i 
[erw exp a (+ gas] ü dH<B,m-ı72 w|exp Ba (1 er me]. SEEN ZEN 
1: \ f 
die Darstellungen | 
Ua | re 
1 on) exp "a 1-+H 3) zen y, BELA NE (29), g 
EM U, y\314 3 Eee Bon 
ol) ex Sarnmel. (Ey e 0), 
1 uoy\ıa De | 
” ze) exp +9 BE). 2. eh, 
1 DV De | 
9 ac \mz exp „> (1+9) np, (BE TR PAD = PN 
Nach Gl. (1),2 G1. (30), und 2. Gl. (32) ist mit einem nach 2Gl. (35) bekannten Fehler die Konstante : 
0 IE, 
p, 7 Bar pP; I ge gs A ee ee NE THE (33). 


3,3 Fehlerabschätzung von 9,, 9, 9, und „Die zugeordnete exakte 5 Sa 
genügt der Differentialgleichung ; s Zu Ya 208 e LOSUNg Pr: 


PniaRt ) (Gr Pr) TIP 0, 


WEENNEEn  Tolimien, Arymptotnche Integration der Störungsdifferentialgleichung 47 
NT ee ee 


während P3 die Ditterentialgleichung 
9%. —IaR(U.— og, —=0.. 
— befriedigt. Also gilt für den Fehler | 
(pn — pa)” —ia R(U —o) (pn —g)" =ia RU) — U) pm... 2.. (89). 
Wie in $2 fassen wir diese Differentialgleichung als ER Differentialgleichung für den 
Fehler auf. Den Fehler wollen wir in einem Intervall abschätzen, das sich von %,, ab ein endliches 


Stück über %. hinaus, etwa bis 2y, erstreckt. In dem Glied auf der rechten Seite von (34) setzen 
wir noch zur Abkürzung ; 


(Dr UN Egg N 435), 


wog beschränkt ist. Da man ein Fundamentalsystem der homogenen Difterentialgleichung in 
1, 9, 9, und o, kennt, so kann man den ns: formal nach der Methode der Variation der Kon- 
stanten hinschreiben: ? 


rs Ä Pe} BY 
—— (Pi — 9) = —P5 ik 9, gY9mdy+ga ) 9, 9Yymdy 
Ua R ; = 2 Ya 2 Ya 
5 ‘ N \ | >(36). 
a SP —e, 9) ypmdy el [99 —9 9 —Y(9P — 9, Ag yııdy 
Ya i Ya \ 
Die unteren Integralgrenzen sind so gewählt, daß der Fehler und seine drei ersten Ableitungen 
gerade bei 2'y, verschwinden. Die Fehlerformel wird auf den Bereich y„ sy s2 4, angewandt, 


ne 
so daß beachtlicherweise „ endlich bleibt. Aus Gl. (10), (26) und (27) folgt 


U— 

Er Keiner De UT 
a RE RT A ER „ [Psı Da—P3Palı ! 2... (37); 

‘ EI 
N AR 1 ka 2 0312 © r) [PsPaa — PP] - - (38). 
Für pr machen wir einen Ansatz ähnlich zu Gl. (27): 
[Ne 2 ; 
Pr = 0-5/4 exp En (1-+) nee TEE ei (BE 


in die Fehlerformel Gl. (36) führen wir die Darstellungen und Ansätze Gl. (10), (26), (27), (29), (30), 
(31), (33),.(37), (38) und (39) ein und benutzen H als Integrationsvariable. Dann ergibt sich: 
H 
RER | 
Du Br 


0 


UN) 
€ Ip | 0- 3/2 Dr 9 Piz Pıı dH 


a 


r 


5 ER 
psp felg a + nme] [a 0222 Voein [—2-2 a +92] menu. 
2H 


—n os Ur "exp 2 1-49) zae| 


3. 
4 (40). 
U’ y\ıua 3, 
| 09/4 e “ 9.eXp ee (1-+ 1) no (Psı Pag — Ps Paı) Pırı AH 
2 & 
ur 5A 5 IR : u. y\ı3a 
— 05/4 eg exp " (1-+:) 432 | F Be RE (P3 Pag — Pr P32) 


2H 


uU y\14 | 
—=n.0—9/4 E (Psı Paa — Ps2 Mi) gexXp = Er 1.4) H9 : Prır A Re 


Das erste der Integrale auf der rechten Seite wird unmittelbar nach der Definition von o_ 3/2 


abgeschätzt, bei den übrigen Integralen wird die zu Gl. (15) analoge Abschätzungsformel für 


H=.2H, 
2H, 


e 


H 


BR f« 


exp | a+nme]an <a, anerp|-a+ nz] . a) 


Für ie FAR Hi: EM a wir zu a. | 6) | 
RER a 19 EL 
NER ne ern 


IN Mag ER RE NL vyıva 
EN ER enge A | Par ya] ausm. 
Be ara Ken on[- Bars m] .. (Me). 
Kan „ Dann erhält man ‚die ee Abschätzungen EB | % Br a U | 
ei SLR AN, N ERW Verla NN BR ee h 
rn N ; Pein.da la Pu a z en 2ER We “ 5 ER 
ä Ve Er Imme—an| lem, Ba NE See (4b), 
RE ; Y = : RN] Nm an | ST alt DREH RE Eh, ge Be 
R wo Be ER rd Pa‘ von & abe Kelnsbante nd er Bar S 
200.2... 83,4 Die Fehlerabschätzung von 9,. 9), p) und m: ; Wachstum von 9 bei y> RE Die. E En 
|  Fehlerabschätzung von 9, und seinen Ableitungen wird RE ähnlich NOTBSROHEBEE: ‚Wir setzen Be: 
N) zu Gl. ur (30), (31) und (32): a Be: 
| Ya u, y \l RS 
9y— en 5) exp |". (1+5) m] Pıy Beh ERROR N: (46a), £ x "A 
hi ff Y 3/4 9 le | Er 
me lg Si er Di me] mv. U Te (ACh), Passage 
1° Papa | Be 
PS (d E- on", (+) Mo& | mes en. 


„(ur (B \ a 
PvE a HA Se exp u (1 g ö). nee | Dil ER. (46.9) 
und lassen die Näherungen mit den exakten Ausdrücken bei Yy übereinstimmen. Bs ergeben sich 7 
die Fehlerschranken: | | 
&. alPalm f BEN.) 
|vıv Pa | Ss Mrpe te Be LEE Eee . (a), 


P : Br 
a Han EN N ER m 


ep RR 
my» —pul < eh ee a7), 


1 |» Im A 
|Bivs Pal ST Tan) RR LEE Ru ans 


diem Gebiet kann man | EUER 
NEE = spatanstapat en 


BL setzen. Deiriid iv: ie man ‚nach a. Sn RT mit pin aus Gl. ) Ai 
“% Je Feanstiman,. muß man Bra Rn 


RE GN u A i RL IE ERREE Mer | 
a r% 27 u 3 ; Fa o m se ” er (Um ga En m + “ I De 
en wählen, wenn man in Er einen relativen Fehler der Größenordnung & zuläßt; denn REN 
er je j gi = io R(Un— - oe Neem a ee “ a ee TE 74 
\ 2 y l J NL ES Dt b So . BA ; r r At Br . Br A Di Betr % Al 2 
n “ F IE i In ca, 22: Ei: Sea | 2 Tre A \ 4 r x - I EH zäh 


-PV enthält demnach einen unbegrenzt lie Bestandtei, wodurch es s für manche Zwecke, 7 
u wird va. 86). 


I 


Bild 5 (links). | 
d?9s /dr 
an’Y dy 
„in Abhängigkeit von H. 


Realteil von 73 


Bild 6 (rechts). 


Ss /ay 
Imaginärteil von dm: AV 


in Abhängigkeit von H. 


$ = 3,5 Veransehaulichung von 2)... Um die ae Funktion @, der Anschauung näher zu 
PER d? 

bringen, sind in den Bildern 5 und 6 der Real- und a von I re 

gegeben, wobei Fı natürlich durch F, angenähert wird. Diese Funktion ist im wesentlichen schon 

von O.Tietgens”’)undH.Schlichting?) tabuliert worden, und zwar besteht der folgende ER 

Zusammenhang zwischen unserer Definition. (links) 2 den Ansätzen dieser beiden Forscher VRR 

| (rechts): ® | Ban... 

x d? dr AT Tee 

e- e 2 za we ee ;)& (mer). Ä N { a 


dn dy B: 


°) Bei der Abfassung dieses As dhnihten verdanke ich Herrn J. Pretsch wesentliche Hilfe. 
’) 0. Tie Eee ns! . Beiträge zum Turbulenzproblem, Diss, Göttingen und Z. angew. Math. Mech. Bd.5 (1925), 
ds - 8.200. 
= ; '®2)H.Schlichting: Amplitndenverteilung und Energiebilanz der kleinen Störungen bei der Plattenströmung. ' 
ey Nachr. d. Ges. d. Wiss, zu Göttingen, Math.-phys. Klasse, Fachgr. I, Bd.1, S.& bis 78, 1935; vgl. besonders $3. 


4 


— Fi wieder- 


sch.ge genau genug. re 1. ‚gerade i in denjenigen 
Verlauf, von @z das, ‚größte, Gewicht ul Beh 
i sr 


IS 3 nz v für die ER), I've I m urn ist, ira u: 2 


a. je ER 2 
RE. ern Son * 
2 SL 4% Are Kuh e . u 


“Fr 


{ ‘ y- Yin + U rar : 
Beis schr kleinen Yu wo. U noch gegen c zu en ist, wird also aus zo 2: 
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TARA RE DE ALPEN 1. Von Ww. Frössel in Göttingen. a 
u . = AR; 


3 In abmeneen mit landleiiz Abkühlu ung‘ treten bei Überschreiten dr ge statt Nebelbil- Bi 
u... dungen Unterkühlungen des Wasserdampfes ein, die erst bei einer bestimmten unter en er stoßartig zusam- Er . 
RN menbrechen und dadurch den Kondensationsstoß erzeugen. : oh 


J 

In exceeding the boundary of saturation under- refrigerations’ of steam instead. of Bas en NESAZ 
are rising in currents of adiabatic refrigeration. These under-refrigerations break down by jerks not IR 
having reached a certain lower boundary and are thus producing the condensing shock. - % FR 


En passant par dessus du point de saturation des sous-refrigerations de vapeur se forment au lieu de con- 
densations brumeuses dans des courants de refrigeration adiabatique. Ces sous-refrigerations s’ecroulent du 
coup seulement a un certain point inferieur et ainsi produisent des coups de conden sation. 2 

B reuennax c anuabeTuyeckuM OXNasKeHHeM, IPH MOCTHKEHHH rpanmısı HPecsimeHng, 
BMeCTO TYMAHHOTO O0Pas0BaHuN, IPOHCXONAT, JIE3KAIMe HUSKe OOBIEHOBEHHOTO, OXIIAKNEHHA 
BOAAHOTO IIAPA, KOTOPEIE HPH MOCTHKeHum OIpeneleHHoH TPAHHIEI TIOHMIKEHNN, TOJIYKO- 
oÖPasHo OCeNalmT, CO3TABAA TAKuM O6Pas0oM 'KOHNEHCAIMOHHBIÄ TOIIYOR. 


nah: 1. Einleitung. Die Messungen von Hochgeschwindigkeitsströmungen werden ‚häufig durch 
ur das Auftreten von Kondensationsstößen !) gestört. Im folgenden wird über derartige Störungen 
Bar“ bei hohen Strömungsgeschwindigkeiten an einer gewölbten Wand berichtet. Die Entstehung 
a der Kondensationsstöße hängt von der atmosphärischen Luftfeuchtigkeit und von den in der R 
Strömung vorkommenden Temperaturen ab. Wie bei der Expansion in Lavaldüsen treten auch. A 
bei der Strömung längs einer Wölbung entsprechend der adiabatischen Abkühlung tiefere Tem- 
peraturen auf, sobald die Schallgeschwindigkeit annähernd erreicht oder überschritten wird. AR SR 
In diesen Gebieten wäre wegen Erreichen der Sättigungsgrenze mit Nebelbildung an Konden- 
sationskernen ?) zu rechnen, Diese kann mit ihrer geringen Nebelbildungsgeschwindigkeit der 
viel schneller stattfindenden Abkühlung in der Strömung jedoch nicht folgen. Statt dessen setzt 
eine Unterkühlung des Wasserdampfes ein, die sich bis zu einer gewissen Grenze fortsetzen läßt, 
um dann bei der starken Übersättigung des Wasserdampfes um so heftiger und stoßartiger zu- 
sammenzubrechen. Diese als Kondensationsstoß bezeichnete lästige Nebenerscheinung bei hohen 
Strömungsgeschwindigkeiten kann in NER ran agL durch Verwendung von trockener Luft 
») Vgl. ron, W, Tollmien $4., 


5 0. Vgl. R. Hermann: Der Kondensationsstoß in Überschall-Windkanaldüsen. Lufo Ba. 19, Lig. 6 (1942), S. 201 
bis 209 


ER 2) Kl. Oswatitsch : Kondensationserscheinungen in Überschalldüsen. Z.angew. Math. Mech. Bd. 22 (1942), 
r.1, 8. 1/14. i 


a a Ein 


be 5 


f Kondensstionsstöße ei kompresibler Srömuns an genölbien 


7? 


\ OnkRmaie bessier, een Im sinne Fall wurde ade Siliea-Gel-Verfahren 3) zur Luft- 
 troeknung benutzt. Eine Verminderung der Luftfeuchtigkeit auf etwa 30% SunEte bereits. zur 
' Vermeidung des Kondensationsstoßes. 

ei 2. Wirkung des Kondensationsstoßes auf die A Zur Untersuchung. 2. Hirt sce: | 
38 von Kondensationsstößen aut die kompressible Strömung längs einer in den Strömungsraum 
VRR vorspringenden Wölbung wurde die im Bild 1 schematisch dargestellte Versuchseinrichtung mit _ 
einem Strömungsquerschnitt von 60 x 60 mm benutzt. Die in Fun eines Kreisbogens vom _ 


5 Radius R=175 mm eingebaute Wölbung war 


rung auf dem Umfang der Wölbung. Die 


 Druckmeßanbohrungen in Abständen von 
20 mm. Sie konnte jedoch selbst um 20 mm 


1 


um ihre Achse drehbar gelagert und ermöglichte 
so beliebige Einstellung der Druckmeßanboh- 


gegenüberliegende gerade Wand besaß feste _ 


stetig verschoben werden, wodurch jede Stelle 


längs der geraden Wand zu erfassen war. Die \  Versuchsstrecke 


Einstellung der Strömungsgeschwindigkeit er- 


h folgte mit dem abflußseitig befindlichen Dit- Bild 1. Schematische Darstellung des Versuchskanals und 


Verstelldiffusors 


fusor. Die Untersuchung der Strömung wurde a Kanalgehäuse mit einem Strömungsquerschnitt 5 


in einer _ Hochgeschwindigkeitsanlage durch- _ von 60 x 60 mm, 


® b drehbar gelagerte Wölbung, Wölbungsradius 
geführt, bei der „atmosphärische Luft in unter Se 295 kim, Köliuneshöher 6 bean 


'brochenem Betrieb in einen luftleer gepump- e Diffusor mit verstellbarem Durchflußquerschnitt. 


ten Behälter strömt. 


; ‘ Die Meßergebnisse werden in den Bildern 2 bis 5 veranschaulicht. In die zwar schema- 
Be aber maßstäbliche, Darstellung des Strömungskanals ist der Kondensationsstoß ge- 


strichelt eingezeichnet. Daran schließt sich der Verdichtungsstoß als dick ausgezogene Linie an. 


' Die Ausbreitung der Stöße senkrecht zur Strömung ist durch Vergleich mit aer kanalhöhe ab- 


schätzbar. Die Stöße selbst sind mit der bekannten Schlierenmethode von 1öpler?) auf 
weißes Papier projiziert und nachgezeichnet. Hier sei bemerkt, daß die Kondensation bei aieser 


‚Schlierenbeleuchtung auch mit dem bloßen Auge sichtbar ist, wenn mian in.Strömungsrichtung 


in die Einsaugdüse hineinsieht. Als bläulicher Dunst hebt sich der Nebel deutlich von der Um- 
gebung ab. Besonders beim Auslaufen eines Versuches ist die stetige Veränderung der Nebelgrenze 
bis zum Verschwinden sehr eindrucksvoll. Gleichzeitig enthalven die Bilder die Auftragung 
der gemessenen Druckverteilungen an der gewölbten, sowie an der geraden Wand, wobeı der 


Ördinatenursprung auf die Mitte der mit aer Kanalwand zusammentallenden Wölbungssehne 


gelegt ist. 3 
Im Bild 2 wird der Strömungsvorgang (Strömungsrichtung von Hanks nach rechts) an der 


| gewölbten Wand mit dem Wölbungsradius £ = 75mm, der Kanalhöhe A = 55 mm und der 


Wölbungshöhe f = 5 mm bei einer Anströmgeschwindigkeit von Ma, = = = 0,68 entsprechend 
v 
einem Druckverhältnis = — 0,732 wiedergegeben. Darin bedeuten Ma, die Machsche Zahl in 


Ik — 
der von der Wölbung ungestörten Anströmung, w, die Strömungsgeschwindigkeit, c, die Schall- 
geschwindigkeit, p, den statischen Druck im gleicnen Raume und 9, den Enwstehungsaruck der 
Strömung (Kesselaruck). In Strömungsrichtung entsteht an der wölbung zuerst eın Konden- 


_ sationsstoß a, der in diesem Bilde von unten nach oben verlaufend bis etwa zur Kanalmitte 


reicht. Daran schließt sich ein zweiter, von oben nach unten fortlaufender Kondensationsstoß b 
an, der an dem folgenden Verdichtungsstoß c endet. Beide Kondensationsstöße werden von der 
gewölbten Wand aurch die Reibungsschicht getrennt, in der infolge der Wandreibung höhere 
‘Temperaturen herrschen und daher aie Voraussetzungen zur Nebelbildung nicht vorhanaen sind. 


Die relative Luftfeuchtigkeit betrug bei den Versuchen in der Atmosphäre etwa 77%. Schon bei 


Verminderung der relativen Luftfeuchtigkeit auf etwa 30% trat der Kondensationsstoß nicht 


“ mehr auf.‘ Die Größe der Ausbreitung quer zur Strömung wie auch die Härte des Konden- 


sationsstoßes ist natürlich weitgehendst von-der atmosphärischen Luftfeuchtigkeit abhängig. 
Die dargestellten Druckverteilungen d und e sind an der gewölbten und an der geraden 
Wand gemessen und aus den unterschiedlichen Bezeiehnungen erkennbar» Während an der ge- 


raden Wand nur ein geringer Druckabfall mit darauf folgendem Druckanstieg als Auswirkung 


der gegenüberliegenden Wölbung stattfindet, erfolgt an der gewölbten Wand anströmseitig.bis 
zum Wölbungsgipfel planmäßig die bekannte Druckabsenkung infolge der an der Wölbung 


3) Die Lufttrocknungsanlage wurde von der Silica-Gel-Gesellschaft, er geliefert. 
“) A. Töpler: Ostwalds Klassiker, Bd.157 (1864). 
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gerade Wand | 


£ 


Bild 4. , Bild 5. 


Bild 2 bis 5. Schematische Darstellung der Strömungsbegrenzung mit eingezeichneten Kondensations- und Verdichtungs- 


stößen und eingetragenen Druckverteilungen. 


a und b Kondensationsstöße, d Druckverteilungen an der gewölbten Wand, 
c Verdichtungsstöße, e Druckverteilungen an der geraden Wand. 


auftretenden Geschwindigkeitserhöhung. Die Stetigkeit des Kurvenverlaufes wird jedoch an der 
Einsatzstelle des Kondensationsstoßes unterbrochen. Hier entsteht eine Verschiebung der Kurve 
nach oben durch Hinzuführen der Kondensationswärme, wie die gemessene Abweichung von der 
ungestört verlaufenden strichpunktierten Kurve veranschaulicht. Der weitere steile Anstieg 
der Druckverteilung ist eine Folge des Verdichtungsstoßes. | 

Der im Bild 3 dargestellte Strömungsvorgang ist durch weitere Öffnung des Diffusors 


(vgl. Bild 1) erreicht worden. Bei geringer Steigerung der Anströmgeschwindigkeit auf 3 


Ma, = 0,698 wird nunmehr unter beibehaltenen Kanalbedingungen der gesamte Strömungsquer- 
schnitt von dem Kondensationsstoß a und 5b durchsetzt, wenn man die warmen Reibungsschichten 
an den Begrenzungswänden unberücksichtigt läßt. Die Wirkung auf die Druckverteilung geht 
aus der Abweichung der Kurve d von der strichpunktiert gezeichneten Fortsetzung des un- 
gestörten Teiles dieser Kurvehervor. Das Störungsgebiet der Kurve wird wieder von der sprung- 
haften Druckerhöhung des stromabwärts folgenden Verdichtungsstoßes c abgeschlossen, der 


ebenfalls bis zur geraden Wand reicht. Die Druckverteilung an der geraden Wand wird in 


diesem Fall auch durch plötzliche Druckanstiege unterbrochen, wovon der erste an der Berüh- 
rungsstelle des Kondensationsstoßes b mit der geraden Wand einsetzt. Darauf folgt wieder ein 
Absinken der Druckkurve bis zum zweiten Druckanstieg, der aber nunmehr von dem Verdich- 
tungsstoß c herrührt. m. | 
Die Bilder 4 und 5 enthalten die Strömungsverhältnisse für die höhere Wölbung f = 10 mm, 
den Wölbungsradiuss R=75mm und die Kanalhöhe k = 60 mm bei den Anströmgeschwin- 


Er; ; > _ ar 
n von geringen Kurvenveränd 
3 im wesentlich 


Bild6b. 
Bild 6a und b. Schlierenaufnahme einer Kondensations- und 
Verdichtungsstoßform mit Schlierenblendenstellung. 
a) parallel zur Strömung, b) senkrecht zur Strömung. 


heller Streifen. Sie wird von dem Kondensationsstoß kaum merklich beeinflußt, wogegen beim 

- Verdichtungsstoß deutlich erkennbar die Ablösung von der Wand stattfindet. An der geraden 
Wand wird die Reibungsschicht entsprechend dem entgegengesetzten Druckgefälle dunkel ab- 
gebildet und verhindert hier die scharfe Abzeichnung der Begrenzungswand. Im Bild 6 b fällt 
die Verschleierung der Wände bei der um 90° gedrehten Schlierenblende fort, da in dieser Stellung 
das Druckgefälle in der Reibungsschicht parallel zur Blende verläuft. Man erkennt auf diesem 

"Bilde auch deutlich, wie voraufgehend erwähnt, daß der Kondensationsstoß wegen der wärme- 
ren Reibungsschichten nicht bis an die Begrenzungswände reicht. 


3. Zusammenfassung. Bei Hochgeschwindigkeitsströmungen mit atmosphärischen An- 
fangszuständen treten unter gewissen Voraussetzungen Kondensationsstöße auf, die unbeachtet 
bei Messungen Anlaß zu Fehlern geben können. Ihre Entstehungsursache ist die hohe Luft- 
feuchtigkeit und die adiabatische Abkühlung in der Strömung, bei der eine Unterkühlung des 
in der Luft befindlichen Wasserdampfes mit darauffolgender stoßartiger Kondensation statt- 
findet. Die dadurch freiwerdende Wärme vergrößert die Gesamtenergie der Strömung, worauf 
die Ermittlung der Strömungsvorgänge aufgebaut ist. Diese störende Erscheinung läßt sich 
durch Verwendung eines Lufttrockners vollkommen beseitigen, 


Eingegangen: Aug. 1944, 
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ebenen P 
veränderlichen 


Taking into consideration the finite tramsverse velocities at the wall appearing in higher degrees of concen- £ 


‚ En ce qui concerne le probleme de la surface plane batgnee d’une manidre longitudinale le champ 


de vitesse et le champ de temperature de valeurs specifiques variables peuvent Etre caleules d’une manibre rela- 


tivement simple au moyen d’une &quation differentielle, demontree par 'E. Pohlhausen et d’une methode 


iterative s’appuyant lü-dessus. Prenant en comsideration les vitesses tramsversales et finies se produisant 


au mur & un. degre de concentration plus elevee, on peut aussi caleuler de cette maniere le champ de diffusion. 2, 


B mpononsHoO LmpoTeramımefi POBHOH INIOCKOCTU TIONe CKOPOCTH H TEeMIEePaTypH 
YCTAHaBIMBaAeTCH IIPH USMEHAWINEHCH OMEHKE BEINecTB CPABHHTENBHO IIPOCTEIM CIIOCOOOM, 
‚ecHu IIPH 9TOM npuMmenäierca yksaHHoe E. IlorpxayseuoMm uHTerpasHoe YpaBHeHme u 0C- 
HOBAHHBIH HA 9TOM ?Ke HHTETPANBHBEIH CIOCOO 0 IpumMeHenHns. TakuMm 3ke IIyTeM BO3MO3KHO 


JCTaHOBUTB Auhhysmonnoe TIONe Ha, CTeHe, IPHHHMAA BO BHHMAHHe HACTymamımue upu 


HOBBIIIEHHBIX KOHNEHTPAUNAX OKOHYATEJIBHLIE TTOTIEPEYHEBIE CKOPOCTH. f 


I. Einleitung. Die laminare Grenzschichtgleichung für die längs angeströmte ebene Platte’ 


bei konstanten Stoffwerten wurde für das Strömungsfeld vonPrandtlfi]undBlasius[2], 


für das Temperaturfeld von E.Pohl-hausen [3] gelöst. Da Temperatur- und Geschwindig- 


keitsfeld zusammenfallen, wenn kinematische Zähigkeit ID und Temperaturleitfähigkeit (a) 


gleich groß sind (Pr = = 1), so ist in der von E. Pohlhausen angegebenen Formel] für das 


Temperaturfeld auch eine Lösung für das Geschwindigkeitsfeld in Form einer Integralgleichung 
enthalten. Pierey und Preston [4] zeigten, daß man, von einer groben Näherung aus- 
gehend, mit Hilfe dieser Integralgleichung und eines Iterationsverfahrens schon in wenigen 
Schritten die bekannte Lösung von Blasius erhält. Dieses Lösungsverfahren besitzt den Vorteil, ein- 
fach zusein undrelativ nur geringen Zeitaufwand zuerfordern. Es ist — wie im folgenden gezeigt 
wird — besonders geeignet für Grenzschichtrechnungen bei veränderlichen Stoffwerten, weileine 
erste und meist recht gute Näherung in der Lösung für konstante Stoffwerte bereits vorliegt. 

Crocco [5] undv. Kärmän und Tsien [6] haben Geschwindigkeits- und Tem- 


_ peraturfeld bei veränderlichen Stoffwerten berechnet. In beiden Arbeiten werden die Differential- 


gleichungen durch Übergang zu neuen Veränderlichen in eine von der sonst bei Grenzschicht- 
rechnungen üblichen abweichende Form gebracht. Crocco erhält zwei simultane Differential- 
gleichungen 2. Ordnung und löst diese für ein Gas mit der Pr-Zahl 0,725 (Luft). v. Kärmän 
und Tsien behandeln den Fall Pr = 1 und haben nur eine Differentialgleichung zu lösen, da 
dann ‘die Temperatur in einfacher Beziehung zur Geschwindigkeit steht. 

Im folgenden wird gezeigt, daß sich eine Reihe von Grenzschichtproblemen an der ebenen 
Platte mit Hilfe dererwähnten Integralgleichungundeines Iterationsverfahrens auf verhältnis- 
mäßig einfache Art lösen lassen, wobei lediglich Quadraturen durchzuführen sind. 

II. Lösung der Grenzschichtgleichung für veränderliche Stoffwerte. Die Grenzschicht- 
gleichungen für Geschwindigkeits- und Temperaturfeld an der ebenen Platte lauten unter Berück- 
sichtigung veränderlicher Stoffwerte mit den üblichen Bezeichnungen [3]: 


%lou) Ale) _ DIR 

2 = De N Se (la) 
or 6) oT 

pour RN, as) Ve : i (2), 


lg 


solange zu. ho ligkeiten 
AS treten, ist ei berechtigt. BE N TREE ee! - 
Ex Bei konstanten Stolfeeiken lassen sich un Sr 2) Aue eine ee a 
u uooslühen, Be 2 Ja den Sr a u ‚und RR N ai 


1 e n Meere een ie ish: der. Gedanke nee daß dieselbe Vereinfachung auch bei ver- & 
RT änderlichen Stoffwerten Bveleh ist.. Wir N ei R KT ERUSE HERRN j 


en eh N ne. 
Be: REN U die Geschwindigkeit. am Rode Yet Grönzechicht, T, und 7, die die Dr 18 
bi E die Temperatur am Rande der Grenzschicht ist. ‚Die Größe Yin der Dimensionslosen € bedeute ! 


. die kinematische Zähigkeit bei der festen Temperatur 7,, für die zweckmäßigerweise die Wand- 


= 52 temperatur (k=0) oder die Temperatur am Rande der Grenzschicht (k= => =) al me a ® 
Ä RER Die ‚Randbedingungen für Strömungs- und Temperaturfeld lauten: 5 et, 
N a DE a Nr oh = 

Be ee” y>0; E00; os el wu 

© . Wir setzen forner, wobei der Index k den ‚Stotiwert ‚bei der Temperatur T, bezeichnet: 
EEE ER eng A | a: = 5 
een eve Fre Ve a 


x Zunächst erhält man aus Fe NR, 
rel: ou BR A ee % ER s 
. damit ergibt sich aus () durch Einführung von er änd (6) in einer a die weitere Rechnung 

it zweckmäßigen Form: . \ 


RR 7) läßt sich folgender: Ma ir @ ping wenn (7) als Differentialgleichung für die ee 


| Größe ge und ferner f vorübergehend als eine bekannte Funktion von & angesehen wird: 


Fr Se Fr 3. 
BF IE)., a I 
ET (00). J(d) = p° de EEE IE A io) 

"Dabei ist über die ee bereits unter Beachtung der Randbedingung (4) ver- 
fügt. Ebenso erhält man für die dimensionslose Bi © den Ausdruck 


—Prn FH dE 
() 


2 RS Bar Bee MR BEN rn ae a RO), 


wobei Pr, die Pr-Zahl mit den Stoffwerten bei > Temperatur T'; ist. 

Für konstante Stoffwerte (= y=x= 1) sind Geschwindigkeits- und Temperaturfeld 
voneinander unabhängig und wir erhalten aus Gl. (9) den bei E. Pohlhausen [3] angegebenen Aus- 
druck für das Temperaturfeld, welcher für Pr= 1 (v = a) die Lösung für das Geschwindigkeits- 
feld darstellt. Während bei bekanntem Geschwindigkeitsfeld die Lösung für das Temperaturfeld 
nach .G].(9) durch einfache Quadratur möglich ist, stößt die Berechnung des Geschwindigkeits- 
feldes zunächst auf die Schwierigkeit, daß in Gl. (8) die noch unbekannte Geschwindigkeit auf der 
rechten Seite in dem Ausdruck für f vorkommt. Bei dem Lösungsverfahren von Piercy und 

- Preston geht man von einer beliebigen Näherung für » aus und berechnet mit dieser L und J(£) 
in Gl. (8); man N so einen verbesserten Wert von w, welcher Ausgangspunkt für einen neuen 


er berechnet 

"Grundlage. ur 

nützt; damit ergab sich die 

RR Fehler in der Schubspannung 
i | , nr Bei veränderlichen Stoffwerte züsch 

Ar Bid. Dieeinzeinen Nüherungen bei Berochn a erläuterten Lösungsschritte. „mathema 

r t teil 'b 

ZIRETR, ET Art“ mit ‚folgenden. Schrißten. ni: 


ann ERetereeN: Kr vereinen: .* EM 
} id Schritt: "Als Ausgangspunkt werden die Deren Lösungen für konstante Stotfwe T 
ke Nu anpeloplinen, sk RAR BLU A 
Kr a) Für das Eecheindgkeitsprdit are Lösung von Blasius 2 nr 


ER > .b) Für das Temperaturfeld die Lösung von E. 'Pohlhausen [3]. % Be 
2. Schrit: a) Berechnung des Geschwindigkeitsprofiles nach Gl. (8), wobei re RR 
RER, N  veränderlichkeit Ber RInENeRDE das Temperaturprofil nach sont ib 
I u - .°... zugrundegelegt ist. 
BE, 'b) Berechnung des a OR nach Gl. (9) mit dem Gr 
profil nach Schritt.2a); EN der Stoffwerte von ‚der‘ Tempera- ar 
Ba: Yo 136 0, bupwie/in\Sehritt 2.) = / a ed}: RL . 
Behr, Das Verfahren wird wiederholt, bis die endgültige Lösung hinreichend genau erreicht ist, ee re 
allgemeinen drei bis vier ‚Schritte genügen. Be. 


Ba . Zunächst: einige allgemeine Ausführungen bezüglich des Einflusses der tobi Joa srieTicn ddr 
0.» Stoffwerte: Die Strömung mit konstanten Stoffwerten kann man als eine erste Näherung ansehen und nach 
den Abweichungen fragen, welche bei veränderlichen Stoffwerten entstehen. Die Güte dieser Näherung hängt re 
offenbar davon ab, welche Temperatur man für die Stoffwerte ‘bei der „isothermen‘“ Strömung annimmt. Wählen ir 
wir die Wandtemperatur oder die Temperatur am Rande der Grenzschicht als Bezugstemperatur für de so- 
therme Strömung, so,ergibt sich auf Grund physikalischer Anschauungen ebenso wie auf Grund der Formeln, _ ’ 
daß innerhalb der Grenzschicht eine Zunahme der Zähigkeit oder der Dichte eine Widerstandsvermehrung 2: EN 

Folge hat; ähnlich bewirkt eine Zunahme der Wärmeleitfähigkeit und der Dichte eine erhöhte Wärmeabgabe. 
Wie groß uber der Einfluß der Veränderlichkeit der einzelnen Stoffwerte ist, hängt vom Verhältnis der Grenz: N Kur. 
schichtdicken des Temperatur- und Geschwindigkeitsfeldes ab?2).- ; 2 


‘ IH. Strömungs- und Temperaturfeld bei zähen Flüssigkeiten. Te Eigenschaften der Stoff- 
Eh zäher Flüssigkeiten entsprechend wurden Geschwindigkeits- und Temperaturfeld unter 
der Annahme berechnet, daß lediglich die Zähigkeit sich mit der Temperatur ändern sollund 
zwar nach folgender Formel: 
7 ir 2: e)‘ 


Ka RE, 
wo b und T, Konstahte sind, welche so gewählt werden, daß die Temperaturabhängigkeit mög- 


lichst gut wiedergegeben wird. Der Index k soll die Werte 0 oder 1 besitzen, je nach der Wahl. 
der Stoffwerte in der Dimensionslosen &. Es wurde b = 3 gewählt (zähes Schmieröl) und die 


D 


SER RE RD: 


beiden Fälle einer beheizten und gekühlten Platte mit u = bzw.8 und Pr,— 12,5 bzw. 100 


1 
berechnet; es handelt sich also um gleich große Temperäturdifferenzen bei derselben Flüssigkeit, 
da ii. jeweils mit den Stoffwerten der Wandtemperatur gebildet ist, Wenn als Bezugstempe- 
&) ratur 7', gewählt wird, erhält man aus Gl. (10) 
| Ar ae ne ER on N EL 
o( we 1)+ 1 
4ı 
Das Ergebnis der Rechnungen nach dem schrittweisen Verfahren des vorhergehenden Ab- 
schnittes ist aus Bild 2 und 3 zu sehen. In beiden Bildern sind die dimensionslosen Wandab- 


ı) Zur endgültigen Lösung des Geschwindigkeitsfeldes in Bild 1 BenDeneIe eine Hilfskraft insgesamt 10 Stunden, 
?) Nach E. Pohlhausen ist das Verhältnis beider proportional —— . 
Ver 


i dabstan Y für beide Abszissen Bisich ist. Nepch der et 
WIZU ung drei‘ Schritte nötig waren, sind noch die isothermen Geschwindigkeits- 5: 

oe ar bzw. k& )ı bei konstanten Stoffwerten bei den Temperaturen T, bzw. T, über den 
 dimensiönslosen Koordinaten eo Bin & a ). .. der Sa ee Bei 


beheizten Platte bei veränderlicher Zähigkeit. / Zähigkeits- 


F matische Zähigkeit bei der Wandtemperatur 7, bzw. der beit IR N \ Rn de KAT Re n 


“heizte Platte und die. Werte 9=14 und 4 0 20 


TEN 
2 \9%x i f 
Bild 2. Geschwiräigkeite.‘ und ee, an einer “Bil a3, Ges van digkeits- un dar emp: a Ei erteilt \ 


EURE Codes Velen (s% 6o) „isotherme‘“ : an der gekühlten Platte. Bezeichnungen nach Bild 2. 
schwindigkeits- und Temperaturverteilungen, », bzw. v, kine- Mut i RUE \ ! 


I STIDSTACNE T, am Rande der renzEehiehi 


{ Hi 


a %o und « die Schubspannung an der Wand und die je Wärmeübergangszahl En —A ‚(0 N) Y ’ 


9% 
(To und (a )o sind die en Werte bei isothermer Strömung mit der Zähigkeit bei der 
Sa ie Lip bzw. Tg. it | Wan 


! OR k “ Tafel 1. 
cn & Pr & = E = 
Ron NE R Klee , 
a i a (Fo)ı (&)o dEyo d&olo 
beheizte WandY . | 0,841.| 120 | 125 | 158:.|. 1,84 
‘ gekühlte Wand 1,08 0,98 | 100 09551, 3 ol 


Es ist leicht einzusehen, daß mit verschwindender thermischer Grenzschichtdicke (Pr so) die 
Schubspannung T, sich von (7,), immer weniger unterscheidet. Die Rechnungen (Tafel 1), be- 
stätigen dies; eine ähnliche einfache Aussage ist bei der Wärmeübergangszahl von vornherein 
nicht möglich; in den gerechneten Fällen ergibt jedoch («), die beste Näherung für a. 


IV. Strömungs- und Temperaturfeld für Pr = 0,7 (Luft) bei Temperaturveränderlichkeit 


‘ sämtlieher Stoffwerte. Im Temperaturbereich —50° bis 140° lassen sich die Stoffwerte der Luft 


durch ıfolgende Formeln darstellen: 


Mi Ki DIENT 0. = R, Me Fu De WE IR 
wobei T die Temperatur in absoluten Graden ist. Mit £ 
Fr 1 
a TE 9 
erhält man- | 
E=9=f1+91— 99% * 
ve 
und ähnliche Ausdrücke für yp und x. 
. Die Rechnungen nach der in Abschnitt II 
angegebenen Methode wurden für eine g0- 


HI Eh 
\ = Shrlar 
Bild 4. Geschwindigkeits- und Temperaturverteilung an 


einer beheizten Platte für Pr = 0,7 (Luft); sämtliche Stoft- 
werte temperaturveränderlich. 


durchgeführt und ergaben im Gchwindigkeile: 
und Temperaturfeld nur mäßige Abweichun- 
gen von der Form bei isothermer Strömung 
(Tafel 2). Zur Untersuchung der: Verhältnisse 
bei höheren Temperaturunterschieden wurde der Fall 7,=20 und T, =,620%°G gerechnet.) Die 
Geschwindigkeits- und Temperaturfelder zeigen nach Bild 4 ne beträchtliche Abweichun- 
gen von der Form bei konstanten Stoffwerten. &, bzw. &, sind mit den Stoffwerten bei den 


3) Bei den „isothermen‘‘ Temperaturprofilen (@), und (©), sind auch die Pr-Zahlen bei den Temperaturen 7’, bzw. 
T, einzusetzen. Z.B, ist in Bild2 Pr, = 12,5 und Ba = 100. 


Ter peraturen 7, bzw. T, 
‘der Grenzschichtdicke; trotzden R 
\ geringe ‚Abweichungen von den Werten 


er 


do 
a 


Ss a: SET REN 0,556 080 er 
re 0 0,420 | 
| r 0,286 0,235 


Nee erklärt sich. daraus! daß Dr Luft die, Zr der Zähigkeit il, E ae il 

einer Widerstandsvermehrung, die Abnahme der Dichte im Sinne einer "Widerstandsvermind 5: 
- rung wirken und sich beide Einflüsse bei Pr = 0,7, wo thermische und Strömunsgrenzschicht-. 
' dicke ungefähr gleich groß sind, praktisch aufheben. Fast ebenso liegen die Verhältnisse beim 
. Temperaturfeld, weil die Wärmeleitfähigkeit eine ähnliche Abhängigkeit von A Teroperaigi 
' besitzt wie die Zähigkeit. - 5% 
In ähnlicher \Weise kann man die Reibungswärme berücksichtigen; in 2) tritt ‚dann auf EN ;% 


| gen: TE & 


‚der rechten. Seite ein Glied «( =) hinzu and man erhält als Lösung: ae 2 Be y 
Be. NN 
a 0. A(E .B Be NER ROT Ss n 
L; R a ( ei ‚ Br & je 2% j “+ B 5 ar 
au-[l Be Sn * a, 
BR, | 5 2 RE AG 2 (12). 
! $ 3 3 > 
REN a. es v2 aba ROG 
SR Dee AN SR MR N N .. 
As se fe U: e \ ar 6 
en | 1 7 rt, 
3 ee je a8 Ana u, ae: 
"TE ; Ä ER 
R Ber Auch hier läßt sich das RR ehe wenn auch die Rechönarböit etwas 
ee. . größer ist. Für,konstante Stoffwerte reduziert sich (12) auf die bereits von E. Eckert) 
A angegebeneLösung. Durch passende Änderung der Integrationskonstanten i in Gl.(12) 1äßt sich 
Bu, auch das Thermometerproblem (verschwindender Temperaturgradient an den Wand) lösen. A 
Ber. - In Hinblick auf die von Crocco für ein Gas mit Pr = 0,725 bereits durchgeführten Rech- # 
REN nungen wurde darauf verzichtet, nach der neuen Methode Zahlenbeispiele durchzurechnen. t 


V. Anwendung auf ein Diffusionsproblem. In derselben Weise wie das Temperaturfeld 
läßt sich beim Diffusionsproblem an der ebenen Platte auch das Konzentrationsfeld berech- 
nen #) 5). Die Differentialgleichung lautet: 

u, Eye Fr a . (3), x 
wobei k die Diffusionszahl und c die Konzentration ist, welche als Menge des Gases bzw. Dampfes Dr 
in der Volumeneinheit definiert ist. Wir wollen hier die Stoffwerte als konstant ansehen, da- 
gegen berücksichtigen, daß bei stärkeren Konzentrationen die Geschwindigkeit v» an der Wand, 
nicht mehr verschwindet, worauf bereits Nusselt [7] hingewiesen hat. Wenn Flüssigkeit von 
einer Wand verdunstet etwa dadurch, daß ein Gas, z. B. Luft, entlang einer benetzten Wand 
strömt, so gelangt ständig Substanz in die Strömung. Man erhält daher an der Wand » (0) >0; 
wenn andererseits Dampf an der Wand kondensiert oder wenn z. B. ammoniakhaltige Luft über 
ein mit Salzsäure getränktes Fließpapier strömt, erhält man » (0) = 0. Die Randbedingung _ 
für v lautet nach den Gleichungen (100) und (101) von [7]: 


I 1 = 
— —|—) —=»() 222.202 nn (AM), 
RR 0) ER 


Po 


4) EB. Eckert hat über eine Lösung dieses Problems auf der Tagung des VDI-Ausschusses für Wartuerorsehune 
1943 in Bayreuth berichtet, wobei ein Näherungsverfahren ähnlich der Methode von K. Pohlhausen für die Strömungs- 
grenzschicht zur Anwendung kam. 


x °) Bine Abschätzung für das vorliegende Problem wurde von G. Damköhler (Z. Elektrochemie [1942], S. 178) 
gegeben. 


=, 


gun 


u 


er wir offenbar nur in Gl. (17) E= 1 zu setzen. 


F 


} Bild 5. Geschwindigkeitsfeld bei Diffusion mit höheren Kon- 
zentrationen, wobei endliche Quergeschwindigkeiten an der 


Wand auftreten [siehe auch Text zu Gl. (14) bis (16)]. 


ation des Gases bzw. Dampfes, für welchen die Wand durchläs 


Ben ano eure ran 
ei 5 30 RR Do 2 a GAR D Be? NER LEN: 
wobei c, und c, die Konzentrationen an der Wand und am Rande der Grenzschicht sind. Wir 
..  „erhaltenähnlich Gl. (5): N HZ RER Br a BR a: 
DENE Ne ee — | wdE+ 1; IM. = — —: Er 00, a) ht 10 
Tr Yozter-joastz]: m Ei 
ee re 2 \ Zi en R Ya: (2-1) Bi g \ 7 
und analog Gl.(9) folgende Lösung für das Konzentrationsfeld: RER 
en Re 
N 0 4 U Re Be 2 
Ba zo=[e a2; ga=2|Hae. EL, Se 
2 I Er ER ko. 0 N \ fi) j RR ; . 
RER TEEDA N; z ae \ 5 ’ = - f 
wobei eine der Pr-Zahl analoge Größe ist. Um die Geschwindigkeit & zu erhalten, brauchen 


Bild 6. "Konzentrationsverteilung zu Bild 5. 


\ 


In dem M steckt der Konzentrationsgradient an der Wand; man kann jedoch Gl. (17) 
zunächst für beliebige M-Werte lösen und mit Hilfe des aus der Lösung erhaltenen Wertes 


für ei die Größe ee 
=) 
: Po 


zentrationsfelder für die gerechneten M- bzw. 
N-Werte sind aus den Bildern 5 und 6, der 
Konzentrationsgradient an .der Wand aus 
Bild7 ersichtlich. M >O bedeutet Verdunstung 
von der Platte, M< 0 Kondensation bzw. 
Absorption an der Platte; für die Größe 
- > wurde der Wert 0,6 gewählt, welcher mit 
guter Annäherung für die Diffusion von Was- 
ser und Ammoniak in Luft gilt ®). Genau 


genommen ist bei den vorausgesetzten hö-" 
 heren Konzentrationen die Dichte und Zähig- 


keit des Gemisches der beiden Stoffe von der 
Konzentration und die Diffusionszahl von der 
Temperatur abhängig. Mit Hilfe der geschil- 
derten Methode können auch derartige Fälle 


ausrechnen. Die Geschwindigkeits- und Kon- 


201 10 


-10L -05 


Bild 7. Konzentrationsgradient an der Wand und Größe M 
in Abhängigkeit von N [siehe auch Text zu Gl.(15) und (16)}. 


°%) Nach Ten Bosch: Die Wärmeübertragung. Berlin 1936, S. 189 und 257. 


wir 


N 
NT 


tigen, wurden angegeben, aber von der Durchführung der dazugehörigen Rechnungen in Hin. _ 

blick auf die bereits vorliegenden Ergebnisse der Arbeit von Crocco abgesehen. Die hier ent- 

- »wickelten Lösungsmethoden wurden schließlich auch auf den Fall der Diffusion von Beimengun- 
gen angewandt, wo bei stärkeren Konzentrationen endliche Quergeschwindigkeiten an der Wand eh 

‚ auftreten. TEN fas 
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KLEINE MITTEILUNGEN 


Über Stabilität von Regelungen mit Nachlauf- 


. zeit. Im folgenden soll als Beitrag zur Stabilitätsfrage 
“ bei Reglerproblemen der Einfluß einer Nachlaufzeit, 


d.h. einer (grundsätzlich stets vorhandenen, oft jedoch 
vernachlässigbar kleinen) Zeitspanne 7, die zwischen 
der Steuerung und der Wirkung derzugeführten Energie 
verstreicht, mit elementaren Methoden untersucht 
werden. 


1. Aufstellung der Gleichung des 
Problems. 


Es sei zur Fixierung der Vorstellung der Einfach- 
heit halber an ein schwingungsfähiges System von 
einem Freiheitsgrad mit geschwindigkeitsproportio- 
naler Dämpfung gedacht (Bild1) auf welches eine 
Kraft P wirkt, die vom Ausschlage x des Systems ab- 


ck 


E- 


hängt, die also durch den Ausschlag & gesteuert wird, 
wobei aber — und das interessiert uns jetzt haupt- 
sächlich — zwischen der Steuerung und Wirkung der 
Kraft eine konstante Verzögerungszeit 7 liegt, d.h. 
die zur Zeit £ auf das System wirkende Kraft Pi) 
hängt von dem Ausschlag %«-r) ab, den .das System 
zur Zeit 6—t hätte; dabei bezeichnen tiefgestellte ein- 
geklammerte Zeiger t, t—t die Zeitaugenblicke. Wir 
machen die paheliegende Annahme einer linearen 
Steuerung (mit p, q als Konstanten) 


Po=rP+ q&t-r). 


Hat das System die Masse m, Federkonstante c und 


Dämpfungskonstante 2x, so lautet die Bewegungs- 
gleichung der Masse 


Mae +2rämn ter =P+g4&-9)*- (l). 


wobei Punkte Ableitungen nach der Zeit hedeuten. 
Mit den Abkürzungen!) | 


1) Hierbei ist c==g angenommen. Auch im Ausnahmefall e =q 
kann man (1) durch einen Ansatz © — f(t) = yin (2) überführen, 
indem man f(t) =&'t (bzw. =&t? oder =&t? in den weiteren 
Ausnahmefällen 2& = — cr und cr? = 2) setzt. Im Falle ce = q 
bleibt, selbst wenn diein Nr. 2 zu besprechenden Lösungen y für & 
beschränkt bleiben, vnicht beschränkt und dasSystemistinstabil. 


3} 


+ do + 2k yo + @! ya) —a yu-n = 0 (2 | 


LER 


Li 


Diese im folgenden näher zu untersuchende Gleichung 
ist eine „‚Funktionaldifferentialgleichung“®); folge —5634 nd 

des zu einem früheren Zeitaugenblick zu nehmenden jschraffierten Bereiche aufsuchen 
 Gliedes ayct-x) ist sie keine gewöhnliche Differen- Nerte di arak ch 
 tialgleichung mehr und hat eine wesentlich größere yeln A mit negativem Realteil besitzt. 
 Mannigfaltigkeit von Lösungen als nur eine zweipara- Diese Bereiche werden von den Bereichen, n 

 metrige Schar, wie sie eine gewöhliche "Differential- Wurzeln A mit positivem Realteil existieren, getrennt 
. gleichung 2. Ordnung aufweist. Man kann den Ver- durch die Grenzkurven, für deren k-a-Werte WurzenA 
lauf von y(i) in einem Intervall der Länge als be- mit dem Realteil 0 vorhanden sind. Für A = 0 erhält 
liebige stetige, am Endpunkt differenzierbare Funk- man die Gerade a = w® (Bild 3). Für A = iv wird de 
_, tion vorgeben und dann y(t) für alle weiteren Zeiten. 
' als Lösung von (2) forsetzeu. \ . Imaginärteil durch. 


at 


B‘ Ä ae @ — a cos» — 0, LER (6) Er 
0. Bei einer gewöhnlichen Differentialgleichung mit 2krtasim-0 } Re 
Bi konstanten Koeffizienten macht man für die Lösung H IR, f ) 
ra .Y(t) einen e-Ansatz. , er" PER SAe Kachı u 94 beschriebertüoGer Kinneh, N ie) 
DLR RL TLE, BER PART { m re ER EEE) N PB a ae NEN 
| DT a er Pe RL ET 
v ; 008% 


2. Diecharakteristische Gleichung. 


und erhält für A ‘eine algebraische Gleichung. Der 
. gleiche Ansatz ergibt jedoch bei der Funktionaldiffe- 

rentialgleichung wegen # 
Yti-ıy= ei u) — cAte— Ar 


für } die transzendente Gleicung 


B+2ki+@®—ae-t=0:: 2: (A). 


kenne reelle Wurzel \- 


Läßt man hier v von — oo bis + oo (oder der Sym- 
. metrie wegen von 0 bis + oo) laufen, so beschreibt die 
zugehörige k-a-Kurve eine unendliche Anzahl von 
| hyperbelähnlichen Ästen, von denen Bild3 für o = T | 


‚die ersten beiden darstellt, und von denen die weiteren 


we N k | Be Bi 


R 
ESS 
& 


rn 


BEER SSIIOSN 


| Br 
| % 
= + X 
F Bild 2. Bild 3. "6 
Diese Gleichung hat 0, 1, 2 oder 3 reelle Wurzeln je Äste die a-Achse immer weiter außen schneiden; diese 
nach den Werten von k, ®®, a, r und im allgemeinen weiteren Äste sind hier ohne Interesse; denn man stellt 
1 unendlich viele komplexe Wurzeln. Für den Spezial- leicht fest, daß von den ‘Gebieten, in die die Grenz- 
{ de, t Y Re kurve die k-a-Ebene zerlegt, nur das in Bild 3 schraf- 
fallı=l,o=- gibt Bild 2 zur Erläuterung ineiner fierte Gebiet frei ist von Wurzeln mit positivem Real- 
‚k-a-Ebene die Bereiche mit 0, 1, 2, 3 reellen Wurzeln, iu a Ba an! El BERN ee 
und zwar ist zugleich in jedem Bereich angegeben, ob M F5 N diese Entechesa ige = @ Y 8 ei | 
“ © die Wurzeln positiv oder negativ sind, z. B. bedeutet + < alte n Pe: had. je as urch. Beirachtung. der Er 
zwei negative und eine positive Wurzel. | Be. Tr Bas: ). 
Besitzt bei fest gewählten Werten von k, w®, a, x die Besonders ‚betont sei, daß das schraffierte Gebiet 
Gl. (4) eine positive Wurzel oder ein komplexes Wur- mit nur Wurzeln von negativem Realteil auch einen 
zelpaar mit positivem Realteil und damit unbegrenzt Zipfel unterhalb der a-Achse , also für negatives k (bei 
anwachsende Partikularlösungen, so ist das zugehörige Anfachung) besitzt. Umgekehrt stellt man fest, daß 
mechanische System instabil. Zur Untersuchung der _ Mur ein ziemlich kleiner sektorähnlicher Teil der oberen 
Wurzeln der „‚charakteristischen, [O1 Be (4) können wir Halbebene kr> 0) zum schraffierten Bereich gehört. 
zunächst 7 = 1 setzen. Das bedeutet keine Beschrän- Be Base a 
® a we uraq = nas man Gen bekannten allohne Verzogerungs- 
kung der Allgemeinheit, sondern lediglich Festlegung zeit; hier bedeutet positives k Dämpfung und Stabilität und nega- 
des Zeitmaßstabes; die Verzögerungszeit wird als Zeit- tives % Anfachung und Instabilität. Für k = 0 und a> 0 zeigt 
B H BEN R ri . eine einfache graphische Darstellung von A? -- w: und act als 
?) Eine ähnliche Gleichung (mit Verzögerungsglied Y(,_ „)in Funktionen von A, daß für «> @: eine positive reelle Wurzel 
der Ableitung) wurde von F. Reinhardt, Der Parallelbetrieb existiert und damit Instabilität besteht. Für k =0 und «<oO 
von Synchrongeneratoren mit- 1 a konstanter ist zwar keine reelleWurzel vorhanden; aber die Gezaden A= const. 
Verzögerungszeit, Wiss. Veröff. emens-Werk, Bd.'18 (1939) . 5 N k Felge 
8.24 bis 44, aufgestellt, aber mit anderen Methoden behandelt, 1 der k-a-Ebene werden für A> + » beliebig flach Zu für 
Im Gegensatz zu dieser Arbeit brauchen hier zur Berechnung der 4 + =) und dringen in jedes der hyperbelähnlichen Gebiete y; 


benötigten Kurven keineWurzeln von transzendenten Gleichungen 
berechnet zu werden. 


der Grenzkurve hinein. 
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Es kann also eintreten, daß gedämpfte Systeme 
(k> 0) bei Berücksichtigung der Nachlaufzeit in- 
stabil werden. 

Man kann innerhalb der schraffierten Gebiete Teil- 
gebiete aufsuchen, in denen älle Wurzeln A einen Real- 
teil kleiner als — u haben, wobei u eine vorgegebene 
Zahl ist. Die Grenzkurve wird dann mit dem Ansatz 
A=— + iv durch die Gleichungen 


2 uv— (u? + @ —v?).tg v 


2 — 
v— utgv . (6) 
„ll Dh) u > 
sinv 
beschrieben. 


3. Direkter Stabilitätsnachweis. _ 
Für Teilgebiete der in Bild 4 schraffierten Bereiche 
läßt sich ein direkter Stabilitätsnachweis erbringen, 
indem man zeigt, daß unter gewissen Bedingungen jede 
Lösung der Gl. (4) beschränkt bleibt für {> oo. Diese 
Teilgebiete sind in Bild 4 eng schraffiert. 
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Weiter sei jetzt |a|< @®< K°. Diese. Ungleichungen 
beschreiben die in Bild 4 engschraffierten ins Unend- 
liche gehenden Halbstreifen. Jetzt ist r > 0. Es sei 
die Zeit >, +r und M eine obere Schranke für 
den Betrag von y(t) im Intervall —r st St, —7, 
also dort |y(t) | SM. Dann kann man in (9) ab- 
schätzen: Pr 


tz 
A Ri } 
y (15) | <& 7 M es (te) — e(k-+r) (6:—8)] ir +E 


Hier ist 
t ! 
JS lea®—n t—9) — er&tn (-9] ds 
0 


1 1 e-(ktnt STR Ar AR 
TEzr karl gan Bere, 
wobei der Rest ; 

j e-k+nt e-k—nt 
(Ele RE RE ET, 


ist, sofern nur t, genügend groß, etwa t> t,ist. Es sei 


Bild 4. 


Die Methode der Variation der Konstanten ergibt?) 
zunächst die allgemeine Lösung der inhomogenen 
Differentialgleichung 


z+2k2z +02 =f(i) 
in der Gestalt 


o..ie.o, Hall ale Nat ie 


t 
2 (t) - [re-a-n e=9—e-@+nü—njdsl (g) 
bo 


+ cjer(k+n) t + Ca e(k—r)i 


wobeir = / k®— w® gesetzt ist, it, ein.beliebiger Zeit- 
augenblick und c,, c, Integrationskonstanten’bedeuten. 
Wir setzen weiterhin i,=0. Mit f(t) =a yt—ı) er- 
hält man für %(t) die mit (2) gleichwertige ‚‚Funk- 
tionalintegralgleichung‘“ 


t 
| 
Ya) = Ir Yes y: 
1) 


| (9) 
h(it)=c,erktn)i+c,erk—nt, 

Nun denken wir uns den Verlauf von y (t) im Inter- 
vall — Tr St sS0 vorgegeben. Dadurch ist eine Funk- 
tion y(t) für alle 2 festgelegt. Wir zeigen, daß unter 
gewissen Zusatzvoraussetzungen % (t) beschränkt bleibt. 

Es sei k > 0, dann haben die Zahlen — (k + r) und 
—(k —r) negative Realteile, undder Bestandteil % (t) 
klingt mit wachsendem i ab. Es gibt kei gegebenem 
&> 0 eine Zeit t,, so daß 


|%(t) |<'e ist für 6> 1. 
Wir betrachten nur noch Zeiten t> t.. 


[Rn t-9) —er(k+n) E-9))ds + %(t) 
mit 


4) Vgl. 2, B. E. Kamke: Differentialgleichungen reeller 
Funktionen, Leipzig 1930, 8.243. Man überzeugtsich am schnell- 
en der Richtigkeit der Lösung (8),indem man sie in (7) ein- 
setzt, 


i, die größere der beiden Zahlen t,, i, und i* eine be- 
liebige Zahl des Intervalles 


u—T = <Su; 
dann folgt 
2r 


w* 
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vet +) Fer 0) 
Hieraus sieht man wegen |@]|< ®® bereits die Be- 
schränktheit von %(t). Im einzelnen kann man dies 
zeigen, indem man sich zunächst & so klein gewählt 
denkt, daß die Zahl 


san dal lalre 
SE ÄST RT 
w“ Zr 


positiv ausfällt. Weiter sei N die größere der beiden 


Zahlen 
€ 


Max 


— TSt<Su—T 


Dann besagt (10) 
2 
weni +.) +2 


n 


und 


yl 


Die Zahl N ist also eine obere Schranke für |y| auch im 


Intervall 4, —r St St, und genau'so auch in jedem 
der folgenden Intervalle 


uKSEISUHTUFTSISHHZT, ..., 
woraus die Beschränkheit von %(t) für alle i folgt. 


Im Falle 0 Sk <wistr reinimaginär,undmankann 
mit derselben Methode mit Hilfe der Abschätzung 


t 
S |esk—r) d—s) — e-(k-+r) (t—s) |ds 
0 


t 
=2f erk t—s)]|sinör(t—s)|ds 
0 


R 
<2fektnds=- HR mit |Rl<e 


a a 


EEE Ve a a 
a AR 


x Ordnung 


het EN von 'y () ) für ut I< k en pn 2 
beweisen. Auch 
der k-a-Ebene sind in Bild 4 eng schraffiert. 


die sich so ergebenden Bereiche in 


‚ Wir fassen die- Ergebnisse zusammen: In den nicht-. 


schraffierten Gebieten in Bild 4 besteht Instabilität, 
“ insbesondere auch für die großen Teilbereiche der - 
. oberen Halbebene k > 0, die ohne Nachlaufzeit stabil 
wären. In den engschratfierten Gebieten besteht Sta- 
 bilität.. In den. weitschraffierten Gebieten ist die Stabi- 
lität durch die vorstehenden Überlegungen noch nicht 


gesichert, sondern nur wahrscheinlich; dort ‚klingen 


. alle eAt-Lösungen ab. i 

Die Betrachtungen lassen. Eh ohne Schwierigkeit 

auf kompliziertere Systeme mit mehreren. Freiheits- 
‚graden übertragen; man erhält dann Funktionaldiffe- 
rentialgleichungen entsprechend höherer Ordnung. 


Hannover. fen . L. Collatz.' 


SHerlelfüng der Leander Hhadrer homogener 
Differentialgleichungen mit konstanten Koeffi- 
.zienten im 
. .charakteristischen Gleichung mittels der Ope- 
. ratorenrechnung. 


Falle mehrfacher Wurzeln der 


Der homogene lineare Differentialausdruck n-ter 


d gel ara 2 
er ey | 


zur a) y(e 


in welchem die a; Konstanten, i die unabhängige Vari- 


able und y eine Funktion von t bedeuten, läßt sich be- 


kanntlich als Produkt von n linearen ee . 


ratoren 1. Ordnung schreiben: 


Ba d d /d 
eek 


Hierbei bedeuten die x; dien Wurzeln der charakteristi- 
schen Gleichung 

an + am— ig m 2er Hay + an=0 (3). 
Die Reihenfolge der Wurzeln x, in 2) ist lol 
die Faktoroperatoren- 


d 
ra) UeE N  eakee, a (4) 


sind also miteinander vertauschbäar. 


Ey =LLL,...-Iny=0: EN 38 (6) 
durch eine Funktion y (t) bereits hinreichend ist, daß 
dieSchlußoperation Nullergibt, daß also: 


Iny = (4) y=0.oder y= On-e’n! (6). 


(O„ eine Integrationskonstante) 
Da wegen der Vertauschbarkeit jeder Faktoropera- 
tor L; in (2) an letzter Stelle TORIEIPLCH kann, so 
ist jeder Ausdruck 
ER ee De EN (7) 
eine partikuläre Lösung von (4). 


a, 


Aus (2) 
“ folgt, daß für das Erfülltsein der Differentialgleichung: 


Sind von ge norsen Be er Re u 


‚gleich, so gibt (7) noch nicht das vollständige Sy-. 


stem der partikulären Lösungen an. Für eine a 


Wurzel Br ist an dann bekanntlich 
te mag, re y. 


ni belie en ig en Koeffizienten C; Lösung der Diffe- - f | 
‚ rentialgleichung. Auch diese, Lösungen, die durch die 


Schreibweise (5) nicht ohne weiteres in Evidenz gesetzt _ 
werden, lassen sich mittels der IPSRITRUEE 
on herleiten. Es ist nämlich: 


li= (ü 


(feine hanige differenzierbare Funktion 72 De \ 
Damit kommt: a 


E 


1dy= u er Re N i 2 
“ en ’ j 4 n A: v 
| REN a a 
zent. em, Guam 10) 
nn .eoore en mt 
d\ } 
RE y 


Liegt nun eine m-fache Wurzel von (3) VOL, so können, 


wir, da es auf die Reihenfolge (der x; in (9) Daahe an- 


kommt, . - 
in-mt+1=m—mt2=. “..=%y (11) 
annehmen. Damit erhalten wir unter Zusammenfassung 
der letzten m Ditierentialoperationen zu en ‚wegen 
nm nen + Di nm +3 n—m + 2% 
d d 
L® = etit, a e (&— Zı)t, a (ts _&)t | 
a a —a,t, ce: 
Bes me ey 


Die Differentialgleichung L(”) y = 0 wird aber bereits 
erfüllt, wenn die Schlußoperation Null liefert: 
dm 
dım ( 
d. h. wenn e” ®n!. y ein Polynom (m—l)-ten Grades 
in 2 mit beliebigen Koeffizienten oder wenn 


| de —=oy 
Yy = (Om De ‚ SH most a (14) 
.. + 0,)* e”® mit beliebigem ( ist. 


Damit ist auch der Fall der Mehrfachwurzeln nfit 
Hilfe der Operatorenrechnung nach dem gleichen Prin- 
zip behandelt wie der der einfachen Wurzeln. 


Dresden, im Juli 1945. W. Buchheim. 


NACHRICHTEN 


Mathematikertagung in Karlsruhe. 
Vom 10. bis 12. April fand eine, insbesöndere den 
Fragen der angewandten Mathematik gewidmete Ta- 
gung deutscher Mathematiker in Karlsruhe statt. Be- 


. reits dort anwesende Herren trafen sich am Mittwoch 


zu einem Begrüßungsabend. Die Tagung war von 
etwa 80 Teilnehmern besucht. Sie begann nach kurzen 
Begrüßungsansprachen eines Herrn der Landesverwal- 
tung und des Herrn Oberbürgermeisters der- Stadt 
Karlsruhe am Donnerstag, früh mit folgenden Vor- 
trägen aus dem Gebiete der angewandten Mathematik: 


Collatz (Hannover), Eigenwertaufgaben bei einer 
Klasse von linearen Integro-Differentialglei- 
chungen. 

Em.de (Pretzfeld), Divergenz und Rotor in nicht- 

. Zlächennormalen Vektorfeldern. 

Wittich (Göttingen), Konforme Abbildung einfach 
zusammenhängender Gebiete. 

Ullrich (Gießen), Über die Abbildung eines Kreis- 
bogenpolygons. 

Meixn er (Aachen), Neuere Ergebnisse über Sphäroid- 
funktionen. 


niet a (x “0. 
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Unger N Zur see Behandlung 
von Anfangswertproblemen bei gewöhnlichen Dif- 
ferentialgleichungen 2. Ordnung. 


Am Nachmittag folgten Vorträge über Rechengeräte 
und aus der Geschichte der ‘Mathematik: 


Dreyer (Darmstadt), Automatisch arbeitende Ge- 
räte für wissenschaftliche Rechnungen. 

Pösch (Weila. Rh.), The Eniac (die große ameri- 
kanische Rechenmaschine). 

Lorey (Frankfurt a. M.), Karlsruher Mathematiker 
vor 1914 und ihre N: achwirkung bis in unsere Zeit. 


Ferner wurde über die künftige Ausbildung der 
Mathematiker an höheren Schulen diskutiert. Für den 
Abend hatte die Stadt zu einer ausgezeichneten Auf- 
führung der Zauberflöte geladen. 

Der Freitag war der Strömungslehre gewidmet; 
er brachte folgende Vorträge: 


Görtler (Oberwolfach), Zur Grenzschichtentstehung 


an Zylindern bei Anfahrt aus der Ruhe, 

Wieghardt (Göttingen), Über einige Unter- 
suchungen an turbulenten Reibungsschichten. 

Fromm (Alzey), ‚Laminare Strömung Newtonscher 
und Maxwellscher Flüssigkeiten. 

Spannhake (Karlsruhe), Gasdynamik ‚in Ro- 
tationshohlräumen mit kleiner Breite. 

Betz (Göttingen), Berechnung von Gasströmungen 
im Bereich der Schallgeschwindigkeit. 

Sauer (Weil a. Rh.), Bemerkungen zur Charakte- 
ristikentheorie der. partiellen  Differentialglei- 

chungen 2. Ordnung. 

Pinl (Oberrahmede), Zur Integrationstheorie adi- 

abatisch-kompressibler Potentialströmungen, 


Schultz-Grunow (Aachen), Theoretische und - 


experimentelle Durchflußbeiwerte bei nicht statio- 
nären Gasströmungen. 

Reutter (Karlsruhe), Über eine angenäherte quasi- 
lineare Potentialgleichung der ebenen kompres- 
siblen Strömung und ihre mittels der Legendre- 
Transformation zu gewinnenden Lösungen. 

Barth (Nonnenborn), Die Bestimmung des Luft- 
widerstandes von schwebenden staubförmigen 
Teilchen. 


Am Nachmittag folgte die allgemeine Mechanik mit 
folgenden Vorträgen: 


Hamel (Landshut), Die Mechanik des allgemeinen 
Kreisels. 

H. Cremer (Aachen), Über den Zusammenhang 
zwischen den Routhschen und den Hurwitzschen 
Stabilitätskriterien. 

H.Cremer (Aachen), Über die durch eine pariodische 
Wechselkraft in einer Platte hervorgerufenen 


"4 Sch en. 
ee, 
L. Cremer (München), Fin einfachen Vorfahnan 205 
DR der Stabilität t linearer Bere 
systeme 
Federn (Darmstadt), Moderne Den des Aus- 
wuchtproblems. _ Fr 
Der Abend vereinigte die Teilnehmer in der Mena, er 
wo sie durch Darbietungen des Hoc er 
(Brahms, Mozart, ) erfreut wurden. 
Der Sonnabendvormittag brachte dann ln Mathe 
Vorträge aus dem Gebiet der angewandten 
matik der Elasto- und Plastomechanik: 


Pickert (Tübingen), Elementare Behandlung hd 
Helmholtzschen Raumproblemes. war 

Merkel (Karlsruhe), Über die Anwendung der affinen. Br 
Abbildung in der Geodäsie. 1% 

Münzner (Göttingen), Über eine indirekte Methode‘ 
zur statistischen Sicherung. 

Ro: 5 rbach (Göttingen), Chiffrierung und Dechiff- 

erung 

We A ner (Heidelberg), Funktionentheoretische Me- 
thode zur Lösung von praktischen Problemen 
der Elastizitätstheorie. 

Kammüller (Karlsruhe), Die Berücksichtigung u. 
von nachträglichen Querschnittsänderungen bei f 
der Berechnung statisch-unbestimmter Systeme 

Karas (Walldorf), Die Ermittlung des Schalen 
druckes von frisch eingefülltem Beton unter Be- 
rücksichtigung der Erhärtung und Reibung. 

Denkhaus (Braunschweig), Die günstigsten Werte 
für das Oberflächendrücken von zylindrischen = 
Drehstabfedern. & 


Am Nachmittag rings noch zwei weitere Vor- 
träge statt: 


Rohrbach (Göttingen), Ergänzungen zum Vor- | 
trag über Chiffrieren und Dechiffrieren. 
Swida (Karlsruhe), Verfahren zur Bestimmung der 
Tragwerkformänderungen im elastisch-plastischen 
Zustand. 


Autorreferate der in Karlsruhe gehaltenen Vorträge 
werden in einem der nächsten Hefte der Z. angew. 
Math. Mech. erscheinen. 


Die Tagung 'gab den auf den “Gebieten der An- 
wendung der Mathematik arbeitenden Mathematikern 
nach langer Zeit wieder die Gelegenheit zu der so sehr 
erwünschten Fühlungnahme und Aussprache. Es ist 
wohlim Sinne aller Teilnehmer, wenn auch hier den 
beiden Karlsruher Mathematikern Herrn Pöschl_ 
und Herrn 'Klotter herzlicher Dank für die Vor- 
bereitung und Durchführung der glänzend gelungenen 
Tagung.ausgesprochen wird. illers. 
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